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A dynamical coupled-channels (DCC) model for neutrino-nucleon reactions in the resonance region 
is developed. Starting from the DCC model that we have previously developed through an analysis 
of nN, 7 N —> 7 tN, r/N, K A, I\ S reaction data for W < 2.1 GeV, we extend the model of the vector 
current to Q 2 < 3.0 (GeV/c ) 2 by analyzing electron-induced reaction data for both proton and 
neutron targets. We derive axial-current matrix elements that are related to the ttN interactions 
of the DCC model through the Partially Conserved Axial Current (PCAC) relation. Consequently, 
the interference pattern between resonant and non-resonant amplitudes is uniquely determined. We 
calculate cross sections for neutrino-induced meson productions, and compare them with available 
data. Our result for the single-pion production reasonably agrees with the data. We also make a 
comparison with the double-pion production data. Our model is the first DCC model that can give 
the double-pion production cross sections in the resonance region. We also make comparison of our 
result with other existing models to reveal an importance of testing the models in the light of PCAC 
and electron reaction data. The DCC model developed here will be a useful input for constructing a 
neutrino-nucleus reaction model and a neutrino event generator for analyses of neutrino experiments. 
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I. INTRODUCTION 


An experimental observation of the leptonic CP violation and a determination of the neutrino mass hierarchy will 
be central issues in forthcoming next-generation neutrino oscillation experiments. Recent findings of relatively large 
$13 [HI have boosted the momentum of the neutrino physics community towards addressing these issues. For a 
success of the next-generation experiments, a more precise interpretation of data will be necessary. This means that 
more precise knowledge of neutrino-nucleus reactions is critically important because neutrinos are detected in the 
experiments through observing remnants of the neutrino-nucleus reactions. 

Neutrino experiments utilize neutrinos in a wide energy range, and therefore the relevant neutrino-nucleus reactions 
have various microscopic reaction mechanisms depending on the kinematics. For a relatively low-energy neutrino 
(E u <> 1 GeV) relevant to, e.g., the T2K [f|, MiniBooNE [i|, and nuPRISM [7[ experiments, the dominant reaction 
mechanisms are the quasi-elastic (QE) knockout of a nucleon, and quasi-free excitation of the A(1232) resonance 
followed bv a decay into a irN final state. For a higher-energy neutrino (2 < E v < 4 GeV) relevant to, e.g., the 
Minerva [&] and future DUNE [9i] experiments, a large portion of data are from higher resonance excitations and deep 
inelastic scattering (DIS). In order to understand the neutrino-nucleus reactions of these different characteristics, 
obviously, it is important to combine different expertise. For example, nuclear theorists and neutrino experimentalists 
recently organized a collaboration at the J-PARC branch of the KEK theory center [l(J 03 to tackle this challenging 
problem. 

In this work, we focus on studying the neutrino reactions in the resonance region where the total hadronic energy W 
extends, mjv + < W < 2 GeV; rriN (m,r) is the nucleon (pion) mass. Furthermore, as a step toward developing a 

neutrino-nucleus reaction model, we will be concerned with the neutrino reaction on a single nucleon. In the resonance 
region, particularly between the A(1232) and DIS regions, we are still in the stage of developing a single nucleon model 
that is a basic ingredient to construct a neutrino-nucleus reaction model. 

First we discuss experimental data that are crucial to determine form factors associated with axial N-N* ( N*: 
nucleon resonance) transitions. Because of small cross sections, neutrino data are rather scarce. Available data are 
from old bubble chamber experiments at the Argonne National Laboratory (ANL) [l2j and the Brookhaven National 
Laboratory (BNL) [lj|; hydrogen and deuterium targets were used in the experiments. The single pion production 
data for E v < 2 GeV are particularly useful, and theoretical models are confronted with the data to fix (or test) 
the strength of the predominant axial 7V-A(1232) transition. However, there has been the well-known discrepancy 
between the two datasets from ANL and BNL by ~10%. This discrepancy is reflected in the uncertainty of the 
axial NA coupling, leading to theoretical uncertainty for neutrino-nucleus reaction cross sections. Regarding this, an 
interesting progress has been reported recently in Ref. [l4| . In the reference, the authors tried to avoid the neutrino 
flux uncertainty of the old bubble chamber experiments. They took advantage of the fact that the ratio cr(l7r)/cr(07r) 
is fairly unaffected by the neutrino flux uncertainty, and that cr(07r) on the deuterium is relatively well understood. 
Here, we denoted the cross section for the single-7r production [no 7r-emission process, mainly QE] by <t(17t) [cr(07r)]. 
Multiplying er(l7r)/<T(07r) from the two experiments by <r(07r) from the GENIE 2.8 [15| . they obtained er(l7r) for the 
two experiments. They found that the newly obtained <t(17t) are both fairly close to the original ANL data. Once this 
result is established, theoretical uncertainty associated with the strength of the axial N A coupling will be significantly 
reduced. 

Another interesting analysis relevant to ANL and BNL data has been conducted by one of the present authors and 
his coworkers in Ref. [16]. They examined effects of the final state interactions (FSI) on cross sections for the single 
pion production off the deuteron. They found that the orthogonality between the deuteron and final pn scattering 
wave functions significantly reduces the cross sections. Thus the ANL and BNL data from deuterium target would 
need more careful analysis with the FSI taken into account. While this kind of reanalysis is important, still the 
available data are rather scarce. It is highly desirable to have new data that are more precise and abundant. There 
is an idea 0 to upgrade the near detector of the T2K experiment to use D 2 O as a target and study this important 
elementary process. 

Regarding theoretical models, several models have been developed for neutrino-nucleon reactions in the resonance 
region; particularly the A(1232) region has been extensively studied because of its importance. Those models can 
be categorized into three classes according to their dynamical contents. The first class gives amplitudes as a sum of 
Breit-Wigner functions that represent resonant contributions. An example of this type developed recently is found 
in Ref. [l8|, and they considered A(1232) 3/2 + , V(1535) 1/2 - , V(1440) 1/2+ and _/V(1520) 3/2“ resonances. The 
second class of models considers tree-level non-resonant mechanisms in addition to resonant ones of the Breit-Wigner 
type. For example, the authors of Refs. ITol l22l derived tree-level non-resonant mechanisms from a chiral Lagrangian, 
and combined them with A(1232) of the Breit-Wigner type. The model of Ref. [20] has been further extended by 
including the V(1520) 3/2“ resonance [23]. Meanwhile, the authors of Ref. [24j developed a model that contains all 
4-star resonances with masses below 1.8 GeV, and included rather phenomenological non-resonant contributions. A 
model of the third class additionally takes account of the hadronic rescattering, so that the unitarity of the amplitude 
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is satisfied. One of the present authors has developed such a model that works for the A(1232) region [2^, (26|. A 
more complete list of existing models for the resonance region can be found, e.g., in Ref. [271 ]. 

Although substantial efforts have been made recently as seen above, there still remain conceptual and/or practical 
problems in the existing models as follows: 

1. We point out that reactions in the resonance region are multi-channel processes in nature. The relevant channels 
such as nN , m tN, t]N, K A, KY are strongly coupled with each other through the strong interaction. Therefore, 
the multi-channel couplings required by the unitarity are essential physics to be considered. However, no existing 
model takes account of this. 

2. The neutrino-induced double pion production over the entire resonance region has not been seriously studied 
previously. As known from the photo- and electron-reactions, cross sections for the double-pion production are 
comparable or even larger than those for the single-pion production around and beyond the second resonance 
region, and a similar tendency is expected for the neutrino reactions. Therefore, it is important to have a good 
estimate for the double-pion production rate. Although some neutrino-induced double-pion production models 
have been developed previously [28i - l30| , they are supposed to work for the threshold region only, a very limited 
kinematics in the whole resonance region. Thus for a practical calculation of neutrino reactions, for example, the 
authors of Ref. [3l| extrapolated a DIS model to the resonance region, and gradually switched on its contribution 
for W > 1.6 GeV, thereby simulating the double-pion contributions. This is not a well-justified procedure, and 
obviously the situation for the double-pion production models is unsatisfactory. 

3. Interference between resonant and non-resonant amplitudes are not well under control for the axial-current in 
most of the previous models. This is due to the fact that the axial-current was not constructed in a manner 
consistent with the ttN interaction model. More detailed discussion on this will be given later in Sec. IIII D 21 

Our goal here is to develop a neutrino-nucleon reaction model in the resonance region by overcoming the problems 
discussed above. In order to do so, the best available option would be to work with a coupled-channels model. In the 
last few years, we have developed a dynamical coupled-channels (DCC) model to analyze nN, jN —> nN, r]N, K A, KY 
reaction data for a study of the baryon spectroscopy [32j . In there, we have shown that the model is successful in 
giving a reasonable fit to a large amount (~ 23,000 data points) of the data. The model also has been shown to give 
a reasonable prediction for the pion-induced double pion productions [33j |. Thus the DCC model seems a promising 
starting point for developing a neutrino-reaction model in the resonance region. At Q 2 = 0, we already have made 
an extension of the DCC model to the neutrino sector by invoking the PCAC (Partially Conserved Axial Current) 
hypothesis [34]. At this particular kinematics, the cross section is given by the divergence of the axial-current amplitude 
that is related to the nN amplitude through the PCAC relation. However, for describing the neutrino reactions in 
the whole kinematical region ( Q 2 ^ 0), a dynamical model for the vector- and axial-currents is needed. 

Here are what we need to do in practice for extending the DCC model to cover the neutrino reactions. Regarding 
the vector current, we already have fixed the amplitude for the proton target at Q 2 =0 in our previous analysis [32| . 
The remaining task is to determine the Q 2 -dependence of the vector couplings, i.e., form factors. This can be done by 
analyzing data for the single pion electroproduction and inclusive electron scattering. A similar analysis also needs 
to be done with the neutron target model. By combining the vector current amplitudes for the proton and neutron 
targets, we can do the isospin separation of the vector current. This is a necessary step before calculating neutrino 
processes. As for the axial-current matrix elements at Q 2 =0, we derive them so that the consistency, required by the 
PCAC relation, with the DCC nN interaction model is maintained. As a result of this derivation, the interference 
pattern between the resonant and non-resonant amplitudes are uniquely fixed within our model; this is an advantage of 
our approach. For the Q 2 -dependence of the axial-current matrix elements, we still inevitably need to employ a simple 
ansatz due to the lack of experimental information. This is a limitation shared by all the existing neutrino-reaction 
models in the resonance region. 

With the vector- and axial-current models as described above, we calculate the neutrino-induced meson productions 
in the resonance region. We compare our numerical results with available data for single-pion, double-pion, and Kaon 
productions. Particularly, comparison with the double-pion production data is made for the first time with the relevant 
resonance contributions and coupled-channels effects taken into account; the previous double-pion production models 
did not include resonant contributions [28|, (2f|, or include a A r (1440) l/2 + resonance contribution only [30j. Also, a 
comparison with other models would be interesting. Thus we will compare structure functions from the DCC model 
with those from the models of Refs. fl8l [35i. [3(|. 

The organization of this paper is as follows: In Sec. El we present cross section formulae for the neutrino-induced 
meson productions. Then in Sec. ED we start with an overview of the DCC model, and discuss the hadronic vector- and 
axial-current amplitudes of the DCC model. We present our analysis of electron-induced reaction data in Sec. lIVI Then 
we present numerical results for the neutrino reactions in Sec. El followed by a conclusion in Sec. IVII Mathematical 
expressions for matrix elements of the non-resonant axial-current, resonant axial-current, and resonant vector-current 
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are collected in Appendices [A] [B] and [Cj respectively. We also tabulate parameters associated with the vector N-N' 
form factors in Appendix [D] 


II. CROSS SECTION FORMULAE 


The weak interaction Lagrangian for charged-current (CC) processes is given by 

£OC = W* J d 3 x(J 00 (x)i0 0„ w + h c , , (1) 

where Gf = 1.16637 x 10 -5 (GeV) -2 and V u d = 0.974 {373 - The leptonic current is denoted by l GG , and is explicitly 
written as 

l° G (x) = Mx) 7/x(l - lb)^ v {x) ■ (2) 

The hadronic current is 

J GC ( x ) = v f (• x ) - K ( x ) > ( 3 ) 

where V+ and A+ are the vector and axial currents. The superscript + denotes the isospin raising operator. For 
neutral-current (NC) processes, the Lagrangian is given by 


£ nc 


Gf 

y/2 


d 3 xJ^ G (x)l NGll (x) , 


( 4 ) 


with the hadronic and the leptonic currents given by, 

J p C (x) = (1 - sin 2 9 w )V 3 (x) - sin 2 6 w V°(x) - A 3 (a;) , (5) 

l ^ C (.x) = - 75)Vv(z) , (6) 

where V)f is the isoscalar current, and sin 2 0w = 0.23 j37j. 

We are concerned with a meson production reaction in neutrino-nucleon scattering v(p u ) + N(pjy) —> l(p{) + f(pf ), 
where the variables in the parentheses are four-momenta of the corresponding particles. A charged lepton (neutrino) is 
denoted by l for CC (NC) reaction. A hadronic final state (/) in our reaction model is one of two-body meson-baryon 
states (ttN,pN, KA, KT,) or three-body nnN states. The double differential cross section with respect to the final 
lepton distribution in the laboratory frame is given with a lepton tensor, L pl/ , and a hadron tensor, , as 

d 3 CTuN —tlf _ GpCy \pi\ J nv w Y,N^f (n\ 

dmsh 47t 2 |p, | ^ ’ v> 


where Y=CCz', CCD, NC for the neutrino CC, antineutrino CC, and NC reactions; Gy = V^ d for Y=CC^, CCD and 
Cy = 1 for Y=NC. We have also used E a as a not ation for the on-shell energy of a particle a. The energy is related 
to the particle mass (m a ) by E a (p a ) = y/m 2 + p 2 . The leptonic tensor is 


LT = pK + PiPt - grip* ■ Pi) ± ie a ^Pu,aPi,p , (8) 

where e 0123 = -fl and +(—) in the last term is for neutrino (anti-neutrino) reactions. The hadron tensor is given by 


s r p f s n 


Here PJ denotes the summation (integral) over the spin (momenta) states of the final hadrons. We used the 

momentum transfer q p defined by q = (uj,q) = p v — pi- The state vector | N) = \N(pn, s z N ,t z N )) denotes the initial 
nucleon state with momentum py and the z-components of the nucleon spin (s z N ) and isospin (t z N ). The state vector 
(f^^\ stands for the scattering final state with the incoming wave boundary condition. Explicit form of the matrix 
element of the hadron current, (/^|./J(0)|iV), will be specified in the next section. 

We express the hadron tensor for a two-body meson-baryon (MB) final state, i/(p,)+-/V(pjv) -A l(pi) + M(k) + B(p), 
in terms of the matrix element of the hadronic current in the center-of-mass frame of the hadronic system (hCM). This 
expression is useful for working with reaction models in which the matrix elements are most easily calculated in the 
hCM. Also, the out-of-plane angle of the meson (4>m] see Fig. [TJ explicitly appears in the cross section formula. We 






5 



FIG. 1. Momentum variables and scattering angles for a neutrino-induced single meson production reaction. 


choose a coordinate system such that the xz -plane coincides with the lepton scattering plane, and the 2 -axis is along 
the momentum transfer q. Regarding the coordinate system in the hCM, the meson is scattered in the xczc ~plane 
and the direction of the zc-axis is chosen along q. The angle between the lepton and hadron reaction planes is denoted 
by <j)M, and is shown along with our coordinate system in Fig. [T] The hadron tensor for a two-body meson-baryon 
final state is given as 

ti^y.a ( 2 tt ) 3 f , n . A a a <7 \k*\E M (k*)E B (p*)E N (p* N ) 

IV - — 2^ J K 

S W> S B 

x (MB^\jJ(0)\N) h cM(MB^\jJmNy h cM ■ (10) 


Here the matrix clement of the hadron current is evaluated in hCM. The initial nucleon state in hCM is | N) = 
| N(px, Spj, tx)), while the final meson-baryon state in hCM, | MB) = \M(k* ,t z M )B(p* ,s z B ,t z B )), has the momentum 
k* ip*), the isospin z-components t z M (t B ), and the spin z-component 0 (s B ) for the meson (baryon). Here we assume 
the meson is a pseudo-scalar particle. The momenta with the asterisk (*) are those in the hCM. Therefore p* = —k* 
and p* N = — q*. The quantity W = y/(q + pn) 2 is the invariant mass of the hadronic system. We have also introduced 
the Lorentz transformation matrix, A J i , that transforms a vector in hCM to that in the laboratory frame. Nonzero 
elements of A J* are explicitly given as 


A/ = A 2 2 = cos 4>m, A] 2 = —A 2 3 


sin (j > M , A 0 ° = A 3 3 


771 tv + OJ 

W 


A 0 — A 3 

yv 3 — 7i o 


M 

W' 


(ii) 


For a two-pion production, i'(p„) + N(pjy) —> l(pi) + ir(pi) + 7r(p 2 ) + N(j>ni), a formal expression of the hadron 
tensor is given as 

j dp 1 dp 2 dp N '6 < ' 4) (q + p N - pi -p 2 - pn') E n ^ PN ' > 

x(7r7rfV(-)|jJ(0)|iV)(7r7riV(-)|jJ(0)|iV)* , (12) 

where |7T7rAf) = |7r(pi, if )tr(p 2 , s z N ,, t z N ,)) is the nnN final state. The isospin states of the pions and the 

spin-isospin states of the final nucleon are denoted as and s z N i,t z N i , respectively. For identical pions, we need a 
symmetry factor, = 1 + St^t^- An explicit formula of the hadron tensor for the two-pion final state in our reaction 
model will be given in the next section. 

Here we also introduce the inclusive cross section in which all the final hadronic states are integrated over. In the 
vanishing lepton mass limit, the inclusive cross section is given in terms of three structure functions, W, Y (W, Q 2 ) 
(*=1,2,3; Y=CO, CCu, NC): 


d 3 CTvN^lX 
dEi dfli 


G 2 f C y 

2n 2 


E?[2W? sin 2 


(h 

2 


Wn COS 2 


0i 


± WY 


E v + Ei 

tun 



(13) 


where Cy has been introduced in Eq. 0 ; the + (—) sign in front of W/ is for neutrino (anti-neutrino) reactions, and 
0i is the lepton scattering angle. The structure functions can be written with the inclusive hadron tensor evaluated 
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in hCM and with Q 2 = — q 2 as follows: 


^i Y = \ (w Y + wl) 


Wl = 


Q 


2\ W 1 


Y 


Q 2 

1/7* I 2 


W 3 y = -^Im[W Y 


xy\ 


w, 


Wlq*» q *-' + 2^Re[Wl) q *» 


and the dimensionless structure functions are conventionally defined by 

F y = m N W Y , F y = ujWJ , F 3 y = 


The inclusive cross section within our DCC model is the sum of the two-body and three-body final states: 

d^&i/N—tlX _ \ " (PVvN^rlMB \ ^ d Cr^N—^lmrN 
dEidfli dEidQi ■“ dEidfli 

MB ttttN 


where the summation of MB (ttttN) runs over all charge states of irN,pN, KA, KT, (7T7riV). 
Finally, we give the differential cross section with respect to W and Q 2 for our later purpose: 

d 2 a vN ^if _ 27 tW d?a uN ^if 

dW dQ 2 2mjv|Pi/||Pz| dEidQi 


(14) 

(15) 

(16) 

(17) 

(18) 


(19) 


III. HADRONIC CURRENT IN DYNAMICAL COUPLED-CHANNELS MODEL 

A. Overview of Dynamical Coupled-Channels Model 

Here we present a brief and overall description of the DCC reaction model [Hj], which is followed by a more concrete 
and practical formulae in the next subsections. We start with an effective Hamiltonian for the meson-baryon system: 

H = Ho + v + T . (20) 

The Fock space of the model consists of meson-baryon states (7rfV, r/N, A'A, A'E and ttttN states) and ’bare’ excited 
states (N*, A, p, a). Here the bare N* state represents a quark core component of the nucleon resonance; it is dressed 
by the meson cloud to form the resonance. We also include 7rA, pN and <jN states as doorway states of the ttttN 
state. The symbol Hq is the free Hamiltonian of the particles, and v is non-resonant interactions among the two-body 
meson-baryon states and 7T7 t states. The non-resonant interactions are based on the s, f, u-channel hadron-exchange 
mechanisms. T represents transitions between bare excited states and two-body states such as A f> ttN. The 
scattering amplitude (T-matrix element) is obtained by solving the Lippmann-Schwinger equation as explained in 
detail in Ref. [38(; the Lippmann-Schwinger equation can be represented diagrammatically as shown in Fig. [2] The 
two-body scattering T-matrix for 7 channels a —> /3 (a,/3 = ttN, pN, A'A, KT,, 7tA, pN, crN) is given by the sum of 
non-resonant and resonant T-matrix as 

(/3|T( ± )(W)|a) = (P\t£l_ ies (W)\a) + (P\4t\W)\a) , (21) 

where the superscript +(—) indicates the outgoing (incoming) boundary condition. Each term of the above equation 
corresponds to the diagram in Fig. [2] in the same ordering. The non-resonant T-matrix is obtained by solving the 
Lippmann-Schwinger equation, 

(/9|^-resWI«) = mV {± \W) + V^\W)G^Uw)t^ n _ ies (W)}\a) , (22) 

where G^ B (W) is Green function of a meson-baryon state; this is the free Green function for the stable channels 
such as MB = ttN, rjN, KA, KT,, while for the unstable MB — ttA, pN,aN channels the Green function contains 
the self-energy to account for the decay into the ttttN channel. The effective two-body interaction is the sum of the 
interaction v and the contribution of particle exchange ’Z-diagram’ that contains the ttttN intermediate state as 

V (±) (W r ) =v + Z^(W). 


(23) 
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M’ 


B' 



M 

B 



Eq. (21) 


AT 


B' 



M 

B 




Eq.(22) 




+ 



( ' * ' ) 


( ' ' ) _1 + 



Eq.(27) 


FIG. 2. (Color online) Diagrammatic representation of the Lippmann-Schwinger equation. The corresponding equation numbers 
in the text are indicated on the right. The symbols T, t, and V in the figures correspond to T^\W) and f r ^j 1 _ res (IE) in 
Eq. (1211) . and V^iW) in Eq. (1231) . respectively. The third line represents the dressed N* decay vertex, (x^|r|2V*), contained 
in Eq. ([261). The double lines represent IV*-propagators and the fourth line corresponds to Eq. ([27[l. 


Introducing a wave operator ' and a scattering state X ( + / ' for outgoing/incoming boundary condition as 

S1 (±) W = 1 + G& B (W)t£l_ ies (W) , (24) 

|xL ±) >=fi (±) |a) , (25) 

the resonant T-matrix is given by 

(M+\W)\a) = ^(x^ ) |r|iV™)^(W)] m ,n(iV:|r| x i + )) , (26) 

m,n 

where the summation is taken over all the considered bare N* and A states labelled by the indices, m and n. The 
N* propagator is denoted by [£)(W)] m>n , and it can have nonzero off-diagonal elements through the rescattering as 

[D-\W)] m , n = {W - m° m )6 m , n + (N^\rrt+Xw)GW B (W)r\N:) , (27) 

for and N* belonging to the same partial wave. In the above equation, is the bare mass for N* r All the 
parameters appearing in the strong interaction Hamiltonian have already been determined in our previous analysis [32j. 

Electroweak meson production reactions on a nucleon (7 + N —> X, e + N —> e 1 + X and v + N —> l + X) within 
the DCC reaction model are described in terms of the matrix element of the hadron current which is, for example, 
the current in Eq. © for the CC neutrino reaction. In the DCC model, the hadron current consists of non-resonant 
meson-production current and bare resonant current as 

m 1 av (p8) 

u Jnon-res ' J res i 

where we have suppressed the label “Y” (=CC, NC) for simplicity, and we will do the same to the symbols and 
in the rest of this section. The currents, jn Qn _ res and j^ es , are the electroweak counterparts of v and T in the strong 
interaction Hamiltonian. Meson-baryon states produced by the current experience multiple rescattering to form final 
states. These scattering processes are described by T matrix elements [Eq. (l2Tj) ] that satisfy the unitarity condition. 
Thus the matrix element of the hadron current between nucleon and meson-baryon scattering state | is given as 

(a^\J*\N) = {xP\jZ oa _ t JN) 

+ T,^ ) \m*J[D(W)Un(K\jL + ^ {+ \W)G^ ) B (W)j^ on _ I jN) , (29) 

m,n 

where the first term on the r.h.s. is purely from non-resonant dynamics while the second term is due to excitations 
of nucleon resonances. This matrix element is the input to calculate the hadron tensor in Eq. ©. A more definite 
expression of the matrix element will be given in the next subsection. We note that not only the on-shell matrix 
elements but also their off-energy-shell behavior are obtained in the DCC model. 


















Finally, a matrix element of the hadron current that contributes to the ttttN final state is given within the DCC 
model used in this work as follows: 

(mrN^lJ^N) = {TrN\T\N)Gi + ^{W){nN { -'>\J ll \N) 

+ (7rN\r\A)Gi + 2(W)(7rA^\J^\N) 

+ (7nT\T\p)G ( ^(W)(pN^\J^\N) 

+ (nn\r\a)G<$(W)(aN^\J»\N) , (30) 

where (MB^\ contains the hadronic rescattering as has been defined in Eqs. (l24| . (|25|) and (|29l) . while (7riV| and 
(7T7r| are non-interacting states. 

Through our previous analysis of the pion- and photon-induced meson production reactions, we have already 
constructed a DCC model for the strong interaction and the electromagnetic current of the proton at Q 2 = 0. Also, all 
the resonance parameters such as the masses and strong decay widths have been extracted from the DCC model j32| . 
In order to extend the model to calculate neutrino-induced reaction cross sections, our main task is to develop an 
axial-current model for j„ on _ res and j^ es . We also need to determine the Q 2 -dependence of the vector form factors 
associated with N-N* transitions by analyzing both electron-proton and electron-neutron reaction data. 


B. Formulae for numerical calculations 

Based on the discussion in the previous subsection, we here present formulae that are practically used in numerical 
calculations. Particularly, we give expressions for matrix elements of the hadron current that are directly plugged 
in the hadron tensor defined in Sec. El In this subsection, all kinematical variables are those in hCM, and thus we 
suppress the asterisk (*) used in the previous section for simplicity. 

Let us consider the matrix element, defined in Eq. (l29l) . for a meson production off the nucleon induced by the 
hadron current J^, Q 2 )\N(—q)), where MB=ttN , rjN , 7tA, pN, aN , KA, ATE. For calculating 

the matrix element, we find it convenience to expand the matrix element in term of the the helicity-L/S'J mixed 
representation; L, S, and J are the orbital angular momentum, total spin, and total angular momentum, respectively, 
for the final MB state. For a detailed discussion on the LSJ- and mixed-representations we employ, see Appendices C 
and D in Ref. [32j]. Thus the matrix element of the hadron current is expressed as follows: 

(M(k; s M s z M t M t z M )B(-k; (<?, Q 2 ) ■ e(A)| N(-q; s N s z N t N t z N )) 

LSJIL* 

x(LL z SS z \Jsx+m(\))YLL*(k)TMB,jN(A,k,q;W,Q 2 ) , (31) 

where q is taken along the z-axis. The spin (isospin) and its ^-component for a particle x are denoted by s x and s z 
(t x and t*), respectively; tj is the isospin of the hadron current. The total isospin of the final MB system is denoted 
by I. The notation (jijf j 2 jf |J 3 J 3 ) stands for the Clebsch-Gordan coefficient. Also, we denoted the polarization 
of the hadron current by A = f, —1,0,1 in the spherical basis (t: time component), and e(A) is the corresponding 
polarization vector; m(A) = 0,—1,0,1 are for A = t, —1,0,1, respectively. We have introduced the matrix element 
in the helicity-L5'J mixed representation denoted by Tmb,jn(X', k, q; W, Q 2 ) in which the label MB is understood to 
include quantum numbers such as L , S, J, and I. The matrix element on the l.h.s. of Eq. (Ellis <MB(-)|J A (0)|JV) hCM 
in Eq. m- Thus we are now left with evaluating the hadron current matrix elements in the hclicity-LS'J mixed 
representation, Tmb,jn, in order to calculate neutrino cross sections using the formulae presented in Sec. El 

Following the manner in the previous subsection, we divide the matrix element into non-resonant and resonant 
parts as 

TMB,jN(A,k,q;W,Q 2 ) = tMB,jN{A 1 k,q;W,Q 2 ) + t^ 1B JN (X;k,q;W,Q 2 ) , (32) 

where is non-resonant amplitude, corresponding to the first term in Eq. E3h and is calculated by 

tMB,jAr(A; k, q; W, Q 2 ) = vmb,jn{ A; k, q\ Q 2 ) 

+ ^ p 2 dptMB,M'B’(k,p- 1 W)GM'B / {p',W)vM'B',jN(X- 1 p,q- 1 Q 2 ), (33) 

M'B' 

where vmb,jn denotes a tree-level non-resonant current matrix element, (MB\j non _ les ■ e(A)|7V) [j n on-res from 
Eq. ES|], projected onto the helicity-LS’J mixed representation. We will specify vmb,jn in the following subsection. 
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Also, is the non-resonant hadronic scattering amplitude expressed with the LSJ representation, and Gm'B' 

is a meson-baryon Green’s function. The second term in Eq. (1321) is the resonant amplitude, corresponding to the 
second term in Eq. (l29l) . and is given by 


tMB,j N ( A; k, g; W, Q 2 ) = ]T f MB ,N^(k; W)[D(W)] mtn T N , tJN (X- q- W, Q 2 ), 

m,n 

f q; W, Q 2 ) = Tjv*,jiv(A; q, Q 2 ) 

+ ^2 p 2 dpr N * t M’B'(p)GM'B’(p,W)tM'B',JN(X- 1 p,q- 1 W,Q 2 ) , 
M’B> J 


(34) 


(35) 


where Tn*,jn denotes a bare iV*-excitation current, (N*\j res -e(X)\N) [j res from Eq. (1251) ]. for which explicit expressions 
are given in Appendix [Bl ( Appendix [Cl) for the axial-current (vector-current). We have also used Tmb.at*, r mb,n* 
and D(W) that are bare, dressed N* —► MB decay amplitudes and a dressed N* propagator, respectively. These 
quantities and (mBjM'B' and Gmb in Eq. (l33l) have been defined and well discussed in Sec. II A of Ref. [32|], and we 
will not repeat it here. 

Before closing this subsection, we derive the hadron tensor for neutrino-induced double-pion productions. Because 
we have specified coupled channels considered in our DCC model, we can manipulate Eq. (1121) to a more definite form 
to be used in actual calculations. As we mentioned, 7rA, pN, oN are the doorway states into the mrN channel in 
the DCC model, and they contribute to the double pion productions as in Eq. (1301) . In addition, we also include a 
mechanism of the hadron current-induced N —> tt + N transition followed by a perturbative N —> nN process for the 
double-pion productions, as in the first term of Eq. (13(71) . In principle, different doorway state contributions shown 
on the r.h.s. of Eq. © can interfere with each other. Also, even within the 7rA doorway state contribution, there is 
still an interference between diagrams in which the final two pions are interchanged. However, we can estimate these 
interference contributions to be small for total cross sections and cPa/dWdQ 2 because of the fact that the Z-diagrams 
in our currently used DCC model give rather small contributions to the ttN scattering [32}. We note however that 
the interference could be more important for differential cross sections with respect to the pion emission angles, which 
we will not calculate in this work, as has been seen in Ref. [H}. Thus we ignore the interference to derive a simpler 
hadron tensor for the double-pion production. For example, the contribution from the 7 tA doorway state to t:\tt^N 
production is given by 

1 M^ bN {t^ b t% b t N t z N \t A tl b +t z N ) 2 T A ^ bN {k]W) 

KbN 2n E A {k) \W - E na {k) - E A {k) - E ATa (fc; W)\ 2 

xWj’ N ^ A , (36) 

where M nb N is the invariant mass of the 7p,7V pair from the A-decay, and k is given by the relation, W = E Va ( k ) + 
M 2 b N + E 7raA (fc; W) is the self-energy of the 7 t 0 A Green’s function as summarized in Appendix A of Ref. [32}; 

r A -> 7 i- t jv(fc; W) = — 2Im[E 7raA (fc; W)}\ +t z N ) 2 T A ^ nbN (k; W) gives the partial width of A into a given 

isospin state of 7r b iV; J]{ab} indicates the sum over the permutations of the two pions, i.e., {ab} = {12}, {21}; the 

symmetry factor r has been defined below Eq. m- The quasi two-body hadron tensor, W / )^’ Ar_>7r “ A , is calculated 
using Eq. m similarly to the other stable two-body channels. A difference from the stable channel case is that k 
determined above is used instead of the on-shell momentum. The hadron tensors for the double-pion productions due 
to the other doorway states can be derived in a similar manner. 


W, 


Y , N TV A —>■ 7T17T 2 N 




= B„ 


•E 

{ab}' 


pW — m-n 




C. Matrix elements of non-resonant currents 

We here specify matrix elements of non-resonant current, Jn 0n _ res in Eq. (I28[) , for meson productions. The matrix 
elements are projected onto the hclicity-LS'J state and plugged into vmb,jn in Eq. (l33l) . As in Eqs. {3|) and ©, the 
current consists of the vector and axial currents. Matrix elements of the non-resonant vector current at Q 2 =0 have 
been fixed through the previous analysis of photon-induced meson-production data, and explicit expressions are given 
in Appendix D of Ref. [32;[. We also need to fix the Q 2 -dependence of the matrix elements to study electron- and 
neutrino-induced reactions, and we use the parametrization given in Appendix A of Ref. [39j . 

Regarding the axial current, we take advantage of the fact that most of our nN —»■ MB potentials are derived from 
a chiral Lagrangian. Thus, we basically follow the way how the axial current is introduced in the chiral Lagrangian: 
an external axial current (ag Xt ) enters into the chiral Lagrangian in combination with the pion field as d^n + /^Oext 
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where f n is the pion decay constant. Therefore, matrix elements of the hadronic axial current are obtained from those 
of 7 tN —> MB by a replacement k M —> */ 7r e A ‘(A), where is the four-momentum of the incoming pion and e M (A) is 
the polarization vector for <ig Xt with a polarization A. We apply this replacement to the nN —► MB potentials of 
the DCC model presented in Appendix C of Ref. [32j. The tree-level axial current matrix elements constructed in 
this way are the non pion-pole part of the axial currents, (MB\A^p tree • e(A)|iV) (i: isospin component), and their 
expressions are presented in Appendix [A] of this paper. By construction, A(^ p tree and the meson-baryon potential v 
satisfy the PCAC relation at Q 2 = — m 2 fiof : 

(MB\q ■ Aj^^JiV) = iU(MB\v\n'N) . (37) 


The axial-current matrix element A 1 ’ 11 in Eqs. © and © is related to the non pion-pole part, A^ p , by 


A*-" 


^NP 


1 


Q 2 + m: 


-q-Al 


NP ! 


(38) 


where the second term is the pion-pole term. In the DCC model of Ref. [32]], we needed to introduce some meson- 
baryon potentials that are not from a chiral Lagrangian in order to fit a large amount of 7r-induced reaction data. 
Although we cannot apply the above replacement to derive the axial currents for those potentials, we can still add 
the corresponding axial currents to maintain the PCAC relation, as given in Eqs. (IADI) . (IA10I) . (IA13I) . and (IA25I) . 

The Q 2 -dependence of the axial-coupling of the nucleon has been studied through data analyses of quasi-elastic 
neutrino scattering and single pion electroproduction near threshold. In the analyses, the axial mass ( Ma ) of the 
dipole form factor Fa(Q 2 ) = 1/(1 + Q 2 /M\) 2 , has been determined to be Ma = 1.026 ± 0.021 GeV (ilj . We employ 
this axial form factor, and assume that all non-resonant axial current amplitudes have the same Q 2 -dependence. 


D. Matrix elements of lV*-excitation currents 

Here, we specify the last piece of the DCC model, Tn*,jn (A; q , Q 2 ) in Eq. (1551) . 


1. Vector current 

The hadronic vector current contributes to the neutrino-induced reactions in the finite Q 2 region. In Ref. (32j, we 
have done a combined analysis of 7riV, yp —> nN , r]N , K A, K E reaction data, and fixed matrix elements of the vector 
current at Q 2 = 0 for the proton target. The bare W./V*-transition matrix elements induced by the vector current 
are parametrized and presented in Appendix [Cj What we need to do is to extend the matrix elements of the vector 
current of Ref. [32j to the finite Q 2 region for application to the neutrino reactions. More concretely, we determine 
Q 2 -dependence of M£[ N ( Q 2 ), Ej^_ N (Q 2 ), S^_ N (Q 2 ), xa 3/2 (Q 2 ) and xa 1/2 (Q 2 ), which we collectively denote by 
F^ n ,(Q 2 ), for the proton and the neutron. (See Appendix [Cl for the definition of the symbols.) This can be done by 
analyzing data for electron-induced reactions on the proton and the neutron, and analysis results will be presented in 
Sec. IIVI Then we separate the vector form factors for N* of / = 1/2 (/: isospin) into isovector and isoscalar parts as 
discussed in Appendix [Cj Regarding N* of / = 3/2 for which only the isovector current contributes, we can determine 
the vector form factors by analyzing the proton-target data. 


2. Axial current 

Because of rather scarce neutrino reaction data, it is difficult to determine N-N* transition matrix elements induced 
by the axial-current. This is in sharp contrast with the situation for the vector form factors that are well determined 
by a large amount of electromagnetic reaction data. Thus, we need to take a different path to fix the axial form factors. 
The conventional practice is to write down a N-N* transition matrix element induced by the axial-current in a general 
form with three or four form factors. Then the PCAC relation, (N*\q- A^ p |./V) = */ 7r (IV*|r|7rW), is invoked to relate 
the presumably most important axial form factor at Q 2 = —m 2 to the corresponding nNN* coupling. The other form 
factors are ignored except for the pion pole term. We then assume A^ p (Q 2 = —m 2 ) ~ A^ p (Q 2 = 0). In the present 
work, we consider the axial currents for bare N* of the spin-parity l/2 ± , 3/2 ± , 5/2 ± and 7/2 ± , and determine their 
axial form factors at Q 2 = 0 using the above procedure. For more detail including explicit formulae, see Appendix iBl 
It is even more difficult to determine the Q 2 -dependence of the axial couplings to N-N* transitions because of the 
limited amount of data. Thus we assume that the Q 2 -dependence of the axial form factors is the same as that used 
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for the non-resonant axial-current amplitudes, i.e., the conventional dipole form factor with Ma= 1-02 GeV. When 
necessary, we can adjust the axial mass for an axial N-N* coupling to fit available data. 

It is worth emphasizing that a great advantage of our approach over the existing models is that relative phases 
between resonant and non-resonant amplitudes are made under control within the DCC model. This is possible in our 
approach by constructing the axial-current amplitudes and nN interactions consistently with the requirement of the 
PCAC relation. As we will see, the resonant and non-resonant contributions to the neutrino reactions are sometimes 
comparable in magnitude and, in such a circumstance, it is essential to have a well-controlled relative phases between 
them to correctly describe the neutrino reactions. We also note that our DCC nN model, on which the axial-current 
is based, has been extensively tested by data in Ref. j32j . 


IV. ANALYSIS OF ELECTRON-INDUCED REACTION DATA 


In this section, we analyze data for electron-induced reactions off the proton and neutron targets to determine the 
Q 2 dependence of F% N ,(Q 2 ), namely Mj± N * (Q 2 ), E{± N ‘ (Q 2 ), S{± N * (Q 2 ), xa 3/2 (Q 2 ), and xa 3/2 {Q 2 ), which are our 
model parameters associated with the vector transition form factors and are defined in Appendix [U] As explained 
in Sec. IIII D 11 for the isospin 7=1/2 nucleon resonances, the analysis of electron-induced reactions on both the 
proton- and neutron-targets is required to decompose the electromagnetic transition form factors into the isovector 
and isoscalar parts. This decomposition is necessary for calculating the CC and NC reactions. Regarding I = 3/2 
nucleon resonances, on the other hand, we determine the F^ N ,(Q 2 ) by analyzing only proton-target reactions data. 
The data we analyze span the kinematical region of W < 2 GeV and Q 2 < 3 (GeV/c ) 2 that is also shared by neutrino 
reactions for E v < 2 GeV. Meanwhile, it is a challenge for the DCC model to predict cross sections of various final 
hadron states in the finite Q 2 region by adjusting only the ’bare’ N-N* transition form factors. Thus the analysis of 
the electron-induced reaction data also serves as a testing ground for the soundness of the DCC model. 


A. Electron-proton reactions 


Among data for electron-proton reactions in the resonance region, those for the single pion electroproductions 
are the most abundant over a wide range of W and Q 2 . Therefore, these are the most useful to determine the Q 2 
dependence of the p-N* transition form factors. The cross sections for p(e, e'n°)p and p(e,e , 7 r + )n have different 
sensitivities to resonances of different isospin state (1/2 or 3/2). The angular distribution of the pion is useful to 
disentangle the spin-parity of the resonances. Based on the one-photon exchange approximation, a standard formula 
of the angular distribution for the single pion electroproduction can be expressed in terms of virtual photon cross 
sections dcrp(Q 2 , W, cos0*)/dfl* (/3 = T, L, LT,TT, LT') for the 7 *N —> nN process in the hCM as, 


d 5 


ep—te'TtN 

dE e idCl e >dQ,t 


+ 'Sf c “ 2<f;+ h • v ' 2e<i ^ e)! Sf sin ‘ 


(39) 


where 


, Q) E e i q 1 

7 “ 2n 2 Q 2 E e 1 - e ’ 


(40) 


and g 7 = (W 2 — m| r )/2mjv, e = [1 + 2(q 2 /Q 2 ) tan 2 (0 e //2)] _1 , the scattering angle of electron 0 e /, the magnitude of 
the virtual photon three momentum g 7 , and the incident (outgoing) electron energy E e ( E e i ) in the laboratory frame; 
h e is the hclicity of the incoming electron; </>* (</>m in Fig. [l]) is the angle between the n-N plane and the plane of the 
incoming and outgoing electrons. The formulae for calculating dap /dfi* from the amplitudes defined by Eq. (1321) are 
given in Ref. [42|. 

The CLAS Collaboration has collected data f~TTil - |77ol | for the single pion electroproduction off the proton in the 
kinematical region of our interest, as shown in Fig. [3] Then they have extracted from the data the virtual photon 
cross sections introduced above [5l|, i.e., dar/dll*. + ed 0 x/dfi*, daLr/d^l*, darr/dfl^, and d<JLT’/d£l*. We fit these 
virtual photon cross sections to determine the Q 2 dependence of the p-N* transition form factors. As seen in Fig. [3] 
the single pion electroproduction data occupy a substantial portion of the relevant kinematical region of W and Q 2 . 
In some kinematical region, however, we still need more data to fix the vector form factors. In particular, data are 
missing for the W j> 1.4 GeV and low-Q 2 region, and the W J> 1.7 GeV and Q 2 <j 2 (GeV/c) 2 . In those kinematical 
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FIG. 3. (Color online) Kinematical region covered by available single pion electroproduction data from the CLAS Collabora¬ 
tion |43H50l |. The red triangle [blue cross] points indicate the kinematical points where data for p(e,e it°)p [p(e, eV + )n] are 
available. At the green square points, both p(e,e'Tv°)p and p(e,e'n + )n data are available. 


region, we fit the inclusive structure functions from an empirical model due to Christy and Bosted [52j . The inclusive 
cross section is defined as 


d &ep->-e'X 

dE e i df2 e / 


r 7 \<J T (W,Q 2 ) + ea L (W : Q 2 )\ 


(41) 


and the proton structure functions Wf m [cf. Eqs. (fbH) and (fl5l) ] are related to the transverse and longitudinal cross 
sections as 


w r = , 

Wr = 4XQi Q +J ‘’ TiW - Q2) + ‘ ,LlW ' Q2)] ■ <42) 


We remark that an analysis of electron reaction data is also interesting in the context of studying the structure 
of nucleon resonances, and we are conducting a fuller and more dedicated analysis of electroproduction data to be 
reported elsewhere [531 ]. 

We have fitted F^ N , ( Q 2 ) to the data at several Q 2 values where the data are available. All the other parameters in 
the DCC model, 406 (115) parameters for hadronic ("f-p) interactions, are fixed as those determined by the combined 
analysis [HJ of nN,"fp —»■ TrN,r]N, KA, KT, data consisting of 22,348 data points. We have successfully tested the 
DCC-based vector current model with the data covering the whole kinematical region relevant to neutrino reactions 
of E v < 2 GeV. Before presenting numerical results of the analysis, we take one more step as follows. In the course of 
the analysis, we have determined the p-N* vector form factors F^ N ,(Q 2 ) at particular Q 2 values where the data are 
available: Q 2 =0, 0.16, 0.20, 0.24, 0.28, 0.30, 0.32, 0.40, 0.50, 0.525, 0.60, 0.65, 0.75, 0.90, 1.15, 1.45, 1.72, 1.80, 2.05, 
2.20, 2.44, 2.60, 2.91, 3.00 (GeV/c) 2 . When the model is applied to calculating neutrino cross sections, we need the 
form factors at arbitrary values of Q 2 . Therefore, it is convenient to parametrize F^ N »(Q 2 ) with a simple polynomial 
function of Q 2 . We approximate the form factors using the following parametrization: 




i^iv.(Q 2 )~E c ^ 2 ) 

n=0 


(43) 


where are constants. We set AC = 5 in Eq. (1431) and require that Cq for the transverse form factors are exactly 
the same as the values determined by the analysis of photo-reactions, i.e., Cq = F^ N ,{Q 2 = 0). Then the other c£, 
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FIG. 4. (Color online) The virtual photon cross section dor/dO ,* + e dai/df!* (/rb/sr) at Q 2 =0.40 (GeV/c) 2 for p(e, e'ir°)p 
(left) and p(e, e'n + )n (ri ght ) from the DCC model. The number in each panel indicates W (MeV). The data are from Ref. [4j] 
for p(e,e'n°)p, and Ref. |45|] for p(e,e'n + )n. 



FIG. 5. (Color online) The virtual photon cross section dor/dfC + eda^/dO,^ (pb/sr) at Q 2 =1.76 (GeV/c) 2 . The data are 
in the range, 1.72 < Q 2 <1.80 (GeV/c) 2 , and are from Ref. [4dl ] for p(e,e'n°)p, and Refs. [4a. [47l for p(e,e , 7r + )n. The other 
features are the same as those in Fig. [4] 


totally 440 parameters, are determined to fit the obtained form factors Fp N ,(Q 2 ). We present numerical values for 
c(( in Table El of Appendix E] We will use the approximate polynomial parametrization of F^ N ,(Q 2 ) in calculations 
presented hereafter. 

Here we present some selected results of the analysis of electron-proton reactions. Among the five differential 
virtual photon cross sections in Eq. (1551) . only dap/dfl^ and dap/dfl^. survive after integrating over the hadronic 
final states. Therefore, they are the most important in view of applications to the neutrino reactions. Thus we 
show a combination of the virtual photon cross sections, dar/dFl*^ + edap/ dfl*, at <3 2 =0.40, 1-76 and 2.95 (GeV/c) 2 
for p(e,e'ir°)p and p(e,e'TT + )n from the DCC model in Figs. HHHl In the same figures, the corresponding data are 
also shown for comparison. The DCC model fits the data reasonably well for both n 0 and n + channels. We also 
show in Fig. [7] our DCC-based calculation of differential cross sections of the inclusive electron-proton scattering 
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FIG. 6. (Color online) The virtual photon cross section dar/dQ^ + edcrn/dn* (pb/sr) at Q 2 = 2.95 (GeV/c) 2 . The data are 
in the range, 2.91 < Q 2 <3.00 (GeV/c) 2 , and are from Ref. [48|] for p(e, e'iv 0 )p and Refs. Phi l47j| for p(e,en + )n. The other 
features are the same as those in Fig. [4] 





FIG. 7. (Color online) Comparison of DCC-based calculation with data for inclusive electron-proton scattering at 
-E e =5.498 GeV. The red solid curves are for inclusive cross sections while the magenta dashed-curves are for contributions 
from the -kN final states. The range of Q 2 and the electron scattering angle ( 6 e /) are indicated in each panel. The data are 
from Ref. IHSl . 


in comparison with data; the single pion electroproduction cross sections from the DCC model are also presented. 
In each of the panels, the range of Q 2 is indicated, and Q 2 monotonically decreases as W increases. The figures 
show a reasonable agreement between our calculation with the data, and also show the increasing importance of the 
multi-pion production processes above the A(1232) resonance region. However, we observe a discrepancy between 
the model and data in W = 1.3 ~ 1.45 GeV in the left panel of Fig. [TJ Because our DCC model gives a reasonable 
fit to the single pion electroproduction data in this kinematical region as shown in Fig. [2 the discrepancy points to 
a problem in our DCC model in describing double-pion electroproduction in this kinematics. By simply adjusting 
the vector form factors, F^ N ,(Q 2 ), we were not able to fit the single-pion and inclusive data at the same time in 
this kinematics. A resolution to the discrepancy in the inclusive cross sections may require a combined analysis 
including double-pion production data. Also, as Q 2 increases, the DCC model starts to underestimate the inclusive 
cross section towards W ~ 2 GeV where the kinematical region is entering the DIS and multi-meson production 
region. Currently available data of neutrino cross sections in the resonance region are not very precise compared to 
the high precision data of the electron-induced reactions. Therefore, for our present purpose of constructing a model 
of neutrino reaction comparable to the current neutrino data, it is sufficient to fit the electron-induced reactions data 
at the level as presented in Figs. 01121 and we refrain from showing the y 2 values for the fits. 
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COS 0/ COS G n * COS d n * COS 0 COS 0 K * COS 0/ COS 8„* COS 0/ COS 0 K * COS 0/ COS 0„* 


FIG. 8 . (Color online) Unpolarized differential cross sections, dcr/dfl* (pb/sr), for 7 n —► 7 r p. The data are from Refs. I58h81| 


B. Photon-neutron and electron-neutron reactions 

Because a free neutron target is not available, “neutron”-target data are extracted from deuteron-target data. 
Effects of the final state interaction and the Fermi motion on the 7 d —> n~pp reaction has been studied in Ref. [55j. 
Here we analyze the data of pion photoproduction on “neutron” available in literature. We analyze unpolarized 
differential cross sections for 777 . —>■ nN from nN threshold to W = 2 GeV, and determine F// n „ ( Q 2 = 0) and the 
cutoffs Ayf [Eqs. (1C 13D - (1C 1 511 ] for 1=1/2 N* states {F// n ,{Q 2 ) = F/// n , (Q 2 ) for 1= 3/2 N* states). A formula to 
calculate differential cross sections from the amplitudes of Eq. (l32l) can be found in Ref. [381, and we will not repeat 
it here. In the finite Q 2 region, we use empirical inclusive structure functions from Ref. |56( as data to determine 
the transition vector form factors F// n ,(Q 2 ). Bosted et al. [57| fitted inclusive electron-deuteron reaction data to 
obtain their model for the inclusive deuteron structure functions, and the inclusive “neutron” structure functions are 
obtained from that by subtracting the proton structure function of Ref. [52| . We use an improved version [56j of this 
“neutron” structure functions. After determining F// n ,(Q 2 ) at Q 2 = 0, 0.20, ..., 3.00 (GeV/c) 2 at every 0.20 (GeV/c) 2 , 
we parametrize them using Eq. (1431) . as we did for the p-N* vector form factors. We present numerical values for c" 
(258 parameters) and those for the cutoffs A y?- (16 parameters) in Tables Hill and HVl of Appendix [D] The following 
results are obtained with this approximate polynomial parametrization. 

We show unpolarized differential cross sections for yn —> nN calculated with the DCC model in comparison with 
data in Figs. [8] and [H and find a reasonable agreement. We also show a comparison of the DCC-based calculation 
with data of differential cross sections per nucleon for the inclusive electron-deuteron scattering in Fig. [10] In this 
calculation, we simply take an average of electron-proton and electron-neutron differential cross sections in the free 
space. As seen in the figures, the calculated resonance peaks are sharper than the data, indicating that the smearing 
due to the Fermi motion and final state interaction needs to be taken into account to obtain a good agreement with 
the data [87]. Finally we note that a more comprehensive analysis including 777 —> nN data is currently underway, 
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FIG. 9. (Color online) Unpolarized differential cross sections da/dfl^ (pb/sr), for 7 n —> n°n. 


The data are from Refs. [UHH]. 



FIG. 10. (Color online) Comparison of DCC-based calculation with data for differential cross sections per nucleon of inclusive 
electron-deuteron scattering at E e = 4.628 GeV. The red solid curves are for inclusive cross sections while the magenta dashed- 
curves are for contributions from the nN final states. The range of Q 2 and the electron scattering angle (9 e i) are indicated in 
each panel. The data are from Ref. [ 86 j. 


which will be reported elsewhere 


V. RESULTS FOR NEUTRINO REACTIONS 


Before discussing neutrino reaction cross sections, first we examine the inclusive structure function F GC defined 
in Eq. dD at Q 2 = 0. At this particular kinematics, the neutrino cross sections are solely determined by F GC . 
F 2 cc (Q 2 = 0) has been often calculated with a PCAC model in which F GC (Q 2 = 0) is related to (experimental) nN 
total cross section by applying the PCAC hypothesis [13, Within our DCC model, we can calculate F GC (Q 2 = 0) 
either with the nN amplitudes or with the hadronic axial-current amplitudes. In Fig. [TTJ we compare F GC (Q 2 = 0) 
calculated by these two ways, and find a good agreement. Some comments are in order. For the CC zz-proton/CC 
z/-neutron process (left panel of Fig. fill) , only 7=3/2 resonances contribute, and the A(1232) dominates F GG in 
low energies. On the other hand, for CC zz-neutron/CC F-proton process (right panel of Fig. fill) , both 7=1/2 
and 3/2 resonances contribute. Thus not only the A(1232) but also 77(1535) 1/2“ and 77(1520) 3/2~ resonances 
in the second resonance region, 7V(1675) 5/2“ and _/V(1680) 5/2 + resonances in the third resonance region create 
characteristic energy dependence of F GG . Although we found the two calculations agree well by construction of the 
model, we still notice that the A(1232) peak from the axial-current amplitude somewhat overshoots that from the nN 
model. Also, there are slight differences in the second and third resonance region in the right panel of Fig. [Tl] These 
differences originate from the fact that the spatial momentum transfer is fixed either at <3 2 =0 or at Q 2 = — m 2 . For 
a more meaningful comparison, we could have compared the two F GG calculated at the same Q 2 . Indeed, we have 
confirmed a significantly better agreement between F GG (Q 2 = —m 2 ) from the axial current amplitudes and those 
from the nN cross sections, which should be the case by definition of the model. Here, we showed F GC that is actually 
used in calculating the neutrino cross sections, and how much it can deviate from those obtained with the nN cross 
sections. Also, we remark that even though most existing neutrino-nucleon reaction models in the resonance region 
have axial-currents based on the PCAC relation, they do not necessarily give F GC (Q 2 = 0) in agreement with those 
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FIG. 11. (Color online) The structure function at Q 2 = 0. calculated with the axial current amplitude (lines) are 

compared with those from n N cross sections via the PC AC relation (points). The upper (lower) line and point are for inclusive 
( 7 tN) final state(s). The left (right) panel is for CC ^-proton/CC F-neutron (CC ^-neutron/CC Fproton) process. 




E v (GeV) 


FIG. 12. (Color online) Total cross sections for the CC v^p (left) and v^n (right) reactions. 


from the nN cross sections; we will come back to this point later in this section. 

Now we present cross sections for the N reactions. With the DCC model, we can predict contributions from 
all the final states included in our model. Also, the DCC model provides all possible differential cross sections for 
each channel. Here, we present total cross sections for the CC N reactions up to E v = 2 GeV in Fig. [T 2 ] we also 
show them in Fig. [13] in log scale to see contributions from all of the final states. For the proton-target, the single 
pion production dominates in the considered energy region. For the neutron-target, the single pion production is still 
the largest, but double-pion production becomes relatively more important towards E v = 2 GeV. The rjN and KY 
production cross sections are 0(1O - 1 -1O -2 ) smaller. 

Next we examine reaction mechanisms of the N scattering. In Fig. 1141 we break down the single-pion production 
cross sections into several contributions each of which contains a set of certain mechanisms. For the proton-target 
process, the contribution from the A(1232) resonance dominates, while the higher N* contribution is very small. 
The A contribution here is the neutrino cross section calculated with the resonant amplitude, Eq. (0), of the P 33 
partial wave only, while the higher N* contribution is from the resonant amplitude including all partial waves other 
than P 33 . The non-resonant cross sections calculated from the non-resonant amplitude of Eq. (l33l) is small for the 
proton-target process. In contrast, the situation is more complex in the neutron-target process where the A gives a 
smaller contribution and both I =1/2 and 3/2 resonances contribute. As can be seen in the right panel of Fig. [lTJ the 
A dominates for E v <, 1 GeV, and higher resonances and non-resonant mechanisms give comparable contributions 
towards E v ~ 2 GeV. This shows an importance of including both resonant and non-resonant contributions with the 
interferences among them under control, as has been stressed in Sec. IIII D 21 Similarly, the contribution of resonant 
and non-resonant amplitudes are shown in Fig. [15] for the two-pion production reaction. Because A(1232) mainly 
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FIG. 13. (Color online) The same as Fig. [12] but in log scale. 




FIG. 14. (Color online) Various mechanisms contributing to v^p —► p ir + p (left) and v^n —► p nN (right). 


contributes below the mrN production threshold and thus gives a small contribution here, the resonant and non¬ 
resonant contributions are more comparable. Still, the figures show that the resonance-excitations are the main 
mechanism for the double-pion production in this energy region. 

Next we compare the CC neutrino-induced single pion production cross sections from the DCC model with available 
data from Refs. [l2|,H3 in Fig. [TO] The left panel shows the total cross sections for v^p —> p~n + p for which A(1232) 
dominates as we have seen in Fig. 1141 If the A(1232)-dominance persists in the neutron-target processes shown in 
the middle and right panels of Fig. [T6l the isospin Clebsch-Gordan coefficients determine the relative strength as 
aiy^n —> p~'K 0 p)/a(v ll p —> p~ir + p) = 2/9 ~ 0.22, and a{y^n —> p~7 r + n)/<j(y^p —> p~TT + p) = 1/9 ~ 0.11. The actual 
ratios from the DCC model are p~-K°p) / a{v^p p~TT + p) = 0.28, 0.27, 0.29, and a{u^n p~TT + n)/a(v f _ l p —> 

p~TT + p) = 0.13, 0.17, 0.21 at ^=0.5, 1, 1.5 GeV, respectively. The deviations from the naive isospin analysis are 
due to the the non-resonant and higher-resonances contributions mostly in the neutron-target processes, as we have 
seen in Fig. [14] The two datasets from BNL and ANL for ryp —> p~n + p shown in the left panel of Fig. QjJ] are not 
consistent as has been well known, and our result is closer to the BNL data [l2j]. For the other channels, our result 
is fairly consistent with both of the BNL and ANL data. It seems that the bare axial 1V-A(1232) coupling constants 
determined by the PCAC relation are too large to reproduce the ANL data. Because axial N-N* coupling constants 
should be better determined by analyzing neutrino-reaction data, it is tempting to multiply the bare axial A-A(1232) 
coupling constants, 5 ajvA(;i 232)> defined in Eq. (IB19I) by 0.8, so that the DCC model better fits the ANL data. The 
resulting cross sections are shown by the dashed curves in Fig. [TgJ We find that —> p~ir + p) is reduced due to the 

dominance of the A(1232) resonance in this channel, while (riy^n —> fi~nN) is only slightly reduced. As mentioned in 
the introduction, the original data of these two experimental data have been reanalyzed recently 0 , and it is pointed 
out that the discrepancy between the two datasets is resolved. The resulting cross sections are closer to the previous 
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FIG. 15. (Color online) Various mechanisms contributing to v^p —► p 7r + 7r°p (left) and v^n —¥ p 7r + 7r p (right). 
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FIG. 16. (Color online) Comparison of the DCC-based calculation (red solid curves) with data for v^p —¥ p~n + p (left), 
—¥ p~n°p (middle) and v^n —¥ p~n + n (right). The DCC calculation with 0.8 x Sana) 1232 ) is also shown (magenta dashed 
curve). ANL (BNL) data are from Ref. fl2 ] (fl3t|). 


ANL data. However, the number of data is still very limited, and a new measurement of neutrino cross sections 
on the hydrogen and deuterium is highly desirable. We also note that the data shown in Fig. [16] were taken from 
experiments using the deuterium target. Thus one should analyze the data considering the nuclear effects such as the 
initial two-nucleon correlation and the final state interactions. Recently, the authors of Ref. |lf| have taken a first 
step towards such an analysis. They developed a model that consists of elementary amplitudes for neutrino-induced 
single pion production off the nucleon [25j, pion-nucleon rescattering amplitudes, and the deuteron and final NN 
scattering wave functions. Although they did not analyze the ANL and BNL data with their model, they examined 
how much the cross sections at certain kinematics can be changed by considering the nuclear effects. They found that 
the cross sections can be reduced as much as 30% for v^d —> p~n + pn due to the NN rescattering. Meanwhile, the 
cross sections for —> p~TT°pp are hardly changed by the final state interaction. It will be important to analyze the 

ANL and BNL data with this kind of model to determine the axial nucleon current, particularly the axial JV-A(1232) 
transition strength. 

Regarding the NC single pion production, we show results in Fig. [T7] In the left panel, we show the cross sections 
for all final charge states. The ratios a(yp —> upn 0 )/a(yp —> vn'K Jr ) ~ cr(yn —> i/mr 0 )/a(i , n —> iypn~) ~ 2 can 
be mostly understood from the isospin Clebsch-Gordan coefficient accompanied by the A — > 7 tN vertex. A slight 
difference between a(up —> ispn 0 ) and a(yn —> unn 0 ) [also between a(vp —> umT + ) and a(yn —> !/p-7r - )] is mostly 
from different interference patterns between the isovector and isoscalar currents. In the middle panel of Fig. 1171 we 
compare the NC 1 m —> vpn~ cross sections from the DCC model with ANL data [90?], and find a fair consistency. In 
closing this paragraph, we compare our DCC-based result for another single meson production, i/n —^ K + A, with 
data in the right panel of Fig. [lT] Although our result undershoots the data, the data are based on statistically very 
limited number of events (3 events) and we still cannot say something conclusive. 

We also compare our calculation of the single pion production with data for differential cross sections with respect 
to Q 2 ( da/dQ 2 ). To make contact with the data in Refs. BiH, we calculate the following flux-averaged cross 
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FIG. 17. (Color online) (left) DCC-based calculation for NC neutrino-induced single pion productions, (middle) Comparison of 
the DCC model and ANL data f9C| for un — upn~. (right) Comparison of the DCC model and BNL data j9l| for un —> pT K + A. 




FIG. 18. (Color online) The flux-averaged (0.5 GeV < E v < 6 GeV) da/dQ 2 for v^p —> p~n + p. (Left) The DCC-based 
calculation (red solid curve) is compared with the ANL data~ fl2il . The DCC calculation with 0.8 x g^///A(i 232 ) is also shown 
(magenta dashed curve). Contributions from high W (>1.4 GeV) are cut in the calculation in order to be consistent with the 
data. (Right) The DCC-based calculation (red solid curve) is compared with the BNL data [Til . 


sections: 


da 

dQ ^ 


dE v 


N ( E v ) da lnnt ] e ] 

^model ( E v ) dQ 2 


(K) 


dE v 


N{E V ) 

(tmod el(-^v) 


(44) 


where N{E V ) is the number of events at neutrino energy E v , and is given in Fig. 6 of Ref. 0 and Fig. 4 of 
Ref. [9l|. The DCC model gives cross sections denoted as a mode \(E u ). In the left panel of Fig. [lHJ we compare the 
DCC-based calculation with the ANL data [12]. We find here again that the DCC model overshoots the ANL data 
for v^p — > p~TT + p. This tempts us to plot da/dQ 2 obtained with the bare axial N- A(1232) coupling constants, 
5ajva(1232)(PCAC), multiplied by 0.8, as we did in Fig. flGH eft). This result is in better agreement with the ANL 
data as seen in Fig. [18] (left). The assumed dipole form factor (Ma = 1-02 GeV) for the bare axial A^-A(1232) vertex 
seems fairl y co nsistent with the data. In the right panel of Fig. [lSl we compare the DCC-based calculation with the 
BNL data [92j]. The Q 2 "dependence of the BNL data with the arbitrary scale is well explained with our DCC model. 

We finally compare our results for double-pion productions with existing data in Fig. 1191 Although there exist 
a few theoretical works on the neutrino-induced double-pion production near threshold [28li30| . our calculation for 
the first time takes account of relevant resonance contributions for this process. The DCC-based prediction is fairly 
consistent with the data in the order of the magnitude. Particularly, the cross sections for v^p —> pi~TT + n°p from the 
DCC model are in agreement with data. However, the DCC prediction underestimates the v^p —> fjrn + n + n data. 
The rather small ratio of a(v^p —> p~Tt + ir + n)/a(v^p p~TT + n°p) ~ 13% at E v =2 GeV from our calculation can be 
understood as follows. Within the present DCC-based calculation, ttttN final states are from decays of the nN and 
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FIG. 19. (Color online) Comparison of the DCC-based calculation with data for v^p —> p 7r + 7r°p (left), v^p —> p n + n + n 
(middle) and v^n —> p~n + n~p (right). ANL (BNL) data are from Ref. M ( 0 )- 
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FIG. 20. (Color online) Contour plots of d 2 a/dWdQ 2 for u^p —> p~n + p (left) and v^n —>- p~nN (right) at E v = 2 GeV. 


of the 7 tA, pN , crN quasi two-body states. For a neutrino CC process on the proton for which hadronic states have 
I = 3/2, the nN, 7rA, pN channels can contribute. Within the current DCC model, we found that the 7 tA channel 
gives a dominant contribution to the double pion productions. Then, retaining only the 7 tA contribution, the ratio 
is given by the isospin Clebsch-Gordan coefficients as, (j{v^p /i _ 7r + 7r + n)/cr(i' A1 p —>■ p~7r + 7r°p) = 2/13 ~ 15%, in 
good agreement with the ratio from the full calculation. With a very limited dataset, we do not further pursue the 
origin of the difference between our calculation and the data. If the double-pion data are further confirmed, then the 
model needs to incorporate some other mechanisms and/or adjust model parameters of the DCC model to explain 
the data. 

Important and characteristic hadronic dynamics changes as W and Q 2 change. Thus, it would be interesting to see 
double-differential cross sections, d 2 a/d,WdQ 2 , as shown in Fig. [521 for the single-pion productions. The prominent 
peak due to A(1232) has a long tail toward higher Q 2 region. For the neutron-target, the resonant behavior in the 
second resonance region is also seen. Similar contour plots are also shown for double-pion productions in Fig. 1211 
Here, the situation is very different from the single pion case, and the main contributors are resonances in the second 
and third resonance regions. 

Since a comparison of the DCC model with other models is interesting, we compare in Figs. [551 and 1551 the structure 
function F^ G {Q 2 = 0) of the DCC model with the model due to Lalakulich et al. jTg} (LPP model), and the Rein- 
Sehgal (RS) model [35], [III. The LPP model consists of four amplitudes of the Breit-Wigner form for A(1232) 3/2 + , 
-/V/1535) 1/2 - , IV(1440) 1/2+ and _/V(1520) 3/2“ resonances with no background. The RS model consists of 18 
Breit-Wigner terms plus a non-interfering non-resonant background of I = 1/2. On the left (right) panel of Fig. 1551 
we show -F 2 cc for the CC reaction on the proton (neutron) going to the ttN final state. From the comparison, a good 
agreement between the DCC model and the LPP model is found only near the A(1232) peak; otherwise they are 
rather different. The RS model rather undershoot the A(1232) peak, as has been also pointed out in Refs. [271 194M96I ]. 
Near the threshold (W ~ 1.1 GeV), meanwhile, F 2 C (Q 2 = 0) of the DCC model is larger than those of the LPP 
and RS models. A similar tendency persists in the inclusive F 2 CC as shown in Fig. [53] For the LPP model, in 

the left panels of Figs. [551 and [551 are the same because only the A(1232) contributes to the proton-target process, 
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FIG. 21. (Color online) Contour plots of d 2 a/dWdQ 2 for v^p — > p n + iv 0 p (left) and v^n — > p n + n p (right) at E v = 2 GeV. 
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FIG. 22. (Color online) Fp c at Q 2 = 0 from contribution of the single pion production only. The DCC model is compared 
with the LPP model due to Lalakulich et al. [l8| and the Rein-Sehgal (RS) model j35|, [3(|. The left (right) panel is for the CC 
U/j, p n ) reaction. 
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and it decays almost exclusively into the ttN state. As discussed earlier in this paper, F GC at Q 2 = 0 is related to 
the 7 tN cross sections through the PCAC relation, and thus is given almost model-independently. We have shown in 
Fig. [TT] that F 2 cg from the DCC axial current model and that from the precise ttN model agree well in accordance 
with the PCAC relation. Therefore, the difference between the DCC model and the LPP and the RS models in 
F^{Q 2 = 0) reveals a consequence of missing the consistency between the axial-current and the nN interaction in 
the latter models. 


VI. CONCLUSION 

In this work, we have developed a dynamical coupled-channels (DCC) model for neutrino-nucleon reactions in the 
resonance region. Our starting point is the DCC model that we have developed through a comprehensive analysis 
of — > jrN,r]N,KA,KY, data for W < 2.1 GeV [32|. The model has also been shown to give a reasonable 

description of ttN —> nnN [33]. In order to extend the DCC model of Ref. [32] to what works for the neutrino 
reactions, we analyzed data for the single pion photoproduction off the neutron, and also data for the electron 
scattering on both proton and neutron targets. Through the analysis, we determined the Q 2 -dependence of the vector 
form factors up to Q 2 < 3 (GeV/c) 2 . By combining the vector form factors for the proton and neutron, we separated 
the vector form factors into the isovector and isoscalar parts; this isospin separation is a necessary step to apply 
the model to the neutrino reactions. We also derived the axial-current matrix elements. An appealing point of our 
approach is that we can derive the axial-current matrix elements that are linked to the ttN potential of the DCC 
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FIG. 23. (Color online) The inclusive structure function Fp c at Q 2 = 0. 


The other features are the same as those in Fig. [22] 


model through the PCAC relation. As a consequence, relative phases between the non-resonant and resonant axial 
current amplitudes are uniquely determined within the DCC model. The <3 2 -dependences of the axial form factors 
are difficult to determine with the available data. Thus we used the same axial form factors for all the axial N-N* 
vertices. Although this prescription is what we can do best for the moment, we hope to improve this in future if more 
data become available. Then, the preparation for calculating the neutrino-induced meson productions off the nucleon 
is completed. 

We have presented cross sections for the neutrino-induced meson productions for E v < 2 GeV. In this energy region, 
the single-pion production gives the largest contribution. Towards E v ~ 2 GeV, the cross section for the double-pion 
production is getting larger to become 1/8 (1/4) of the single-pion production cross section for the proton (neutron) 
target. Because our DCC model has been determined by analyzing the ttN,jN —> ttN, r]N, KA, KT, data, we can 
also make a quantitative prediction for the neutrino cross sections for pN , K A, and KT, productions. We found that 
cross sections for pN, K A and K E productions are 10 - 2 -10 -3 times smaller than those for the single pion production. 
We have compared our numerical results with the available experimental data. For the single-pion production, our 
result, for which the axial N-N* couplings are fixed by the PCAC relation, is consistent with the BNL data for 
v^p —> 7 T + p, while fair agreement with both ANL and BNL data is found for the neutron target data. Through the 

comparison with the single pion production data for W < 1.4 GeV for which the A(1232)-excitation is the dominant 
mechanism, we were able to study the strength and the Q 2 -dependence 0 f the axial AT-A(1232) coupling. We also 
calculated double-pion production cross sections by taking account of relevant resonance contributions for the first 
time, and compared them with the data. We found a good agreement for v^p —> fi~n + TT°p and u^n —> /j,~tt + tt~p^ but 
not for v^p —> fx~TT + n + n. Because the data are statistically rather poor, it is difficult to make a conclusive judgement 
on the DCC model. We hope new high quality data will become available in near future. We examined reaction 
mechanisms contributing to the single-pion production. For the proton target where only I = 3/2 states contribute, 
the A(1232) dominates. However, for the neutron target where both I = 1/2 and 3/2 states contribute, the A(1232), 
higher N* resonances, and non-resonant mechanisms give comparable contributions towards E v ~ 2 GeV. Thus it 
is very important to have interference patterns among those different mechanisms under control. In this regard, our 
DCC approach has an advantage over the other existing models, as mentioned in the above paragraph. 

Finally, we make some remarks on our future prospect. Our DCC model should be smoothly connected to a DIS 
model in the resonance-DIS overlapping region. Because the DCC model still has degrees of freedom to vary the 
Q 2 -dependence of the axial form factors, we can adjust them to fit the W- and Q 2 -dependences of inclusive cross 
sections from the DIS model in the overlapping region. This needs be done in a future work. Having developed 
the DCC neutrino-nucleon reaction model that covers the whole resonance region, next task should be developing a 
neutrino-nucleus reaction model in which the DCC model describes the elementary processes. The simplest case is 
the neutrino-deuteron reactions, and for that, we can do a relatively well-established quantum mechanical calculation, 
as done in Ref. [n|. It will be interesting to extend the work of Ref. [jjJ by replacing the elementary amplitudes 
therein with those from the DCC model developed here. Then, ANL and BNL data should be reanalyzed with 
this model that takes care of not only the Fermi motion but also the final state interactions (FSI). The importance 
of FSI has been demonstrated in Ref. [l6j]. Also, this development can serve as a preparation for a possible T2K 
experiment 0 that utilizes heavy-water (D 2 O) as the target. Application to heavier nuclei will also be very important. 
In the A(1232) region, the well-developed A-hole model [93, [9g may give a hint to address pion productions in the 
neutrino-nucleus reaction. For even higher energy neutrino reactions, a fully quantum mechanical calculation seems 
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formidable. Combining the elementary amplitudes of the DCC m odel with a hadron transport model may be a 
possible and practical option, as has been done in the literature inm. 
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Appendix A: Matrix elements of non-resonant axial currents 


We explained in Sec. IIII Cl how to derive matrix elements of non-resonant axial-currents to be implemented in the 
DCC model. Here, we present tree-level, non pion-pole (NP) part of the matrix elements for a single meson production, 
{M(k', j)B(p')\A l (q)\N(p)}, where i and j are isospin components and the other variables in the parentheses are four- 
momenta carried by the corresponding nucleon (TV), baryon (B), meson (M), and axial-current (A^); labels for the 
spin and isospin states for the baryons are suppressed. The following expressions in Appendices lAllCl are those in the 
center-of-mass frame of the hadronic system (hCM). The matrix elements can be expressed in the following form: 


{M{k'j),B{ V ')\A^ PMee {q) ■ e(A)| N(p)) 

= (SF X \j E„(- q) Zt -£) ' 


(Al) 


where e At (A) is the polarization vector for the hadron axial current with a polarization A. The Dirac spinor for 
the baryon B is denoted by ub(p) that is also supposed to implicitly contain the isospin spinor. In the following 
subsections, we present expressions for A(n) for each process labeled by the index n. A(n) is composed of several 
terms, denoted as A™, Aetc. for which we also present diagrammatic representations in TABLE [U It is noted 
that, in evaluating the time component of four-momenta contained in the prop agato rs in the following equations, 
we follow the definite procedures defined by the unitary transformation method [10lL Il02j ; see also Appendix C of 
Ref. [H]. 

For the axial-current matrix elements shown in the following subsections, we include at each vertex a dipole form 
factor in the same way as those used for ttN —> MB potentials in Ref. [32|. For a meson-baryon-baryon vertex, we 
include a form factor of the form 

< A2 > 

where k and A are the meson momentum and the cutoff, respectively. For a N-B transition vertex induced by the 
axial current, we also include a form factor of Eq. (IA2I) in which k is chosen to satisfy 


W = \Jk 2 + ml + \Jk 2 + m 2 N , (A3) 

in order to satisfy the PC AC relation, Eq. (1371) . In a t-channcl diagram, there is a meson-meson transition vertex 
induced by the axial-current, and we use a dipole form factor of Eq. (IA2I) where k being the momentum of the 
exchanged meson. For a contact term, we use double dipole form factor, F(k', A)F(fc, A), where k' is the outgoing 
meson momentum. In addition to this hadronic form factor, we also include the axial form factor of the dipole form, 
as discussed in Sec. IIII Cl to take care of the Q 2 -dependence of the axial couplings. 

In the following expressions, we use notations such as k = p — p' and k that satisfies Eq. (IA3I) . We also denote for 
e M (A) for simplicity. The *-th component of Pauli matrix that acts on the nucleon isospin spinors is denoted by t\ We 
denote the pseudoscalar, vector, and scalar mesons by P , V, and S , respectively, and the octet and decuplet baryons 
by B. With this notation, the PBB ', VBB', and SBB' coupling constants are denoted respectively by fpBB 1 , 9vbb' > 
and gsBB while the PP’V and PP'S coupling constants are denoted by gpp'v and gpp'S, respectively. The tensor 
coupling constant of a VBB' coupling is denoted by Ky- For numerical values of couplings, masses, cutoffs appearing 
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TABLE I. Diagrammatic representations for non-resonant axial-current matrix elements. A", A £,... correspond to expressions 
presented in Eqs. (IA5l) - (IA45ll . The wavy lines represent the external axial-currents. 



below in this appendix, we use those determined in Ref. [32j, and listed in TABLEs XI-XIII of the reference. In 
addition, we use f n = 93 MeV. 
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where S^ AI3 is the positive energy part of the A propagator; in the frame where the A is at rest, 

s£ )mn (p) = j; 0 1 + 7 ° (3 5 mn - a m a n ) . (A12) 

o p u — rriA 

Also, the operator T ( T t) generates an isospin transition from I = 1/2 to 3/2 (/ = 3/2 to 1/2) states. We note here 
that we include the cross diagram of s-channel A(1232)-resonance diagram as a part of the non-resonant mechanism, 
as in Eq. (1A71) . For S 31 partial wave, wg also add 
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where e*Jf is the polarization vector for the p-meson, and 


rp ' ~ T 


Sp'* + ^-Mp* ft'~ ft' ip'*) 


N 


and e 0123 = +1 convention is taken. 
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For the K production amplitudes in this and the following subsections, the isospin operator r acts on 
of the initial nucleon and the final K. 
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where the suffix j is the isospin state of E. 
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Appendix B: Bare N* axial current matrix elements 


We present the axial N-N* transition matrix element for each of l/2 ± , 3/2 ± , 5/2 =1= , 7/2 ± bare N* in the helicity 
basis. The matrix element is parametrized in terms of form factors that can be determined at Q 2 = 0 by invoking 
the PCAC relation to the 7r + N —> N* matrix element. Before presenting the expressions for the axial N-N* matrix 
elements, we give the definition for the axial N-N* matrix elements and how they are connected to the 7r + IV —> N* 
matrix elements through the PCAC relation. We also specify how the axial N-N* matrix elements presented in this 
section are implemented in the formulae in Sec. IIII Al 

A tree-level s-channel bare N* amplitude for an axial-current induced single pion production is given in the plane 
wave basis by 
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where the normalization of the amplitude is the same as that of Eq. dm The axial N-N* and N* —»■ nN matrix 
elements are (N*\A l N ,(q) ■ e\N(p)) and {n(k', j)N(p')\hisr*->nN\N*), respectively. Let us first present expressions for 
the N* —y nN matrix element that is subsequently related to the axial matrix element by the PCAC relation. The 
N* —> nN matrix element (n(k', j)N(p')\liN*^nN\N*) is parametrized by 
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= —ig*N,N*{k'){\j-t z N \T N >T^ [ ,)'^2,{LL z -s z N \JN>MN*)YLL^{k') , (B2) 


where (ljl/2f^ r |Tjv*T^ r ») is the isospin Clebsch-Gordan coefficient. The N* has the spin Jn* and the parity (—l) i+1 , 
and it decays into the nN state with the orbital angular momentum L. The vertex function is denoted by g*N,N* ( k') 
that is related to in Eq. (® by 


5ttJV,JV* (*') = i 


. (2n) 3 E N (-k')2E n (k') 


m N 


Ttt N,N*(k') 


I E N (-k')2E n (k') 


L N 


a 


77iV(L,S=l/2),jV* ^ A 2^ + fe /2 


A 2 

ly N* 


(2+L/2) n k '\\ L 


(B3) 


where we have used the parametrization for T w n,n* defined in Ref. [32]. We use numerical values for the coupling 
CkN(l,s=i/ 2 ),n* an d the cutoff Ajv* presented in Ref. [32j |. 

Now we discuss the axial N-N* transition matrix element. First we separate the dependence on the isospin 
components from the matrix element by 
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where is the spin-parity of N*. The non pion-pole part of the matrix element A^ ± (g) can be determined by the 

ttNN* coupling constant g-nN,N* using the PCAC relation, (iV*|g • \N) = if v {N*\hitN^N*\n l N). Since the 

7 tNN* coupling constants have been determined at Q 2 = —m 2 from our analysis of the nN reaction data [H}, we 
determine A(/ p N , (g) at Q 2 = — m 2 using the PCAC relation: 
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for J N . = L + 1/2, and 
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for J N » = L — 1/2; ± is for s z N = ±1/2. As will be shown in the following subsections, we use Eqs. (1B6I) and (1B7I) to 
fix form factors g jN *(Q 2 = — m 2 ) contained in A/) ± (g). Then we take the PCAC hypothesis: g jN *(Q 2 = —to 2 ) ~ 
± ± ^ N * 

g jN *(Q 2 = 0). For the Q 2 -dependence of g J ”*{Q 2 ), as discussed in Sec. IIIID21 we assume the dipole form factor 
that is implicit in the expressions in the following subsections. We note that the axial form factor g J *r* ( Q 2 ) fixed by 
Eqs. (IB6I) and (IB7I) at a range of q acquires a W-dependence that is the same as that of g^ nn* (q); q and W are 
related by W = \fq 2 + ml ± sfq 2 ± m 2 N . Thus we denote the axial form factors by g J N * (W, Q 2 ). The axial-current 
amplitudes determined in this way satisfy the PCAC relation with the nN amplitudes at any W, as shown in Fig, fill 
Finally, the axial-current part of the vertex function r^r* jn in Eq. (1351) is related to the axial-current matrix element 

by 
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1. Spin 1/2 N* 

The matrix element of axial vector current between the nucleon and spin l/2 ± N* [Aj± (g) • e in Eq. M is 
generally given, the induced tensor term being omitted, by 
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where the upper (lower) operator is for positive (negative) parity AC*; g 1 / 2 ( Q 2 ) and gp 2 ( Q 2 ) are form factors. The 
matrix element for the divergence of the axial current is 
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where we have introduced the IF-dependence in the form factors. The PCAC hypothesis dictates g]( 2 ( W,Q 2 = 
—m 2 ) ~ gT (W, Q 2 = 0). The gp term can be understood as the pion pole term, and thus is given by 
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where the helicity of the axial current is indicated by e^) with A = +1,0, t in the spherical basis. The sign + is for 
= ±1/2. For the 1/2“TV* state, the helicity amplitudes are 
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2. Spin 3/2 IV* 


The matrix element of axial vector current between the nucleon and 3/2 ± N* \Aj± ^(q) in Eq. m] is generally 
given by 
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where we have introduced the bF-dependence in the form factors. The PCAC hypothesis dictates g^ 2 ( W,Q 2 = 
—ml) ~ g g / 2 (W. Q 2 =0). From the PCAC and the pion dominance, we have a pion-pole term: 
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3. Spin 5/2 A* 

The matrix element of axial vector current between the nucleon and 5/2 ± N* \Aj± ^(q) in Eq. (EDI is generally 
given by 
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where we have introduced the IF-dependence in the form factors. The PCAC hypothesis dictates g\j 2 (W,Q 2 = 
—m 2 ) ~ g^ 2 (W. Q 2 =0). From the PCAC and the pion dominance, we have a pion-pole term: 


gl /2± (W,Q 2 ) = • 

Considering g 3 and 174 terms, the helicity amplitudes for 5/2+TV* are 
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4. Spin 7/2 A* 

The matrix element of axial vector current between the nucleon and 7/2 ± TV* \Aj± ^{q) in Eq. M is generally 
given by 
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where gy/ 2 (Q 2 ) (n = 1...4) are form factors. The matrix element of the divergence of the axial current is 
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where we have introduced the IF-dependence in the form factors. The PCAC hypothesis dictates gj 2 (W,Q 2 = 
— m 2 ) ~ g^ 2 (W. Q 2 =0). From the PCAC and the pion dominance, we have a pion-pole term: 
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Considering g% and g 4 terms, the helicity amplitudes for 7/2 + are 
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Appendix C: Bare A* vector current matrix elements 


Following the formulation of Ref. 
representation as 


we parametrize a bare 7 (*^N —> N* vertex matrix element in the helicity 
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where A 7 and Ajv are the helicity quantum numbers of the (virtual) photon and the nucleon, respectively, and q R and 
qo are defined by Mjv* = q R + E^{q R ) and W = qo + E^{—q), respectively. For N* of I = 1/2, this vertex can be 
separated into the isovector (F JV ) and isoscalar (r /s ) parts as follows: 
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where (10l^|l^) is the isospin Clebsch-Gordan coefficient. Regarding N* of I = 3/2, they have only the isovector 
part as given by 
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The quantities introduced above (rjy) ^jy, rj)f. yAr ) correspond to the vector part of Tn*,jn in Eq. (l35l) . The helicity 
amplitudes G% N (q) in Eq. (IC1I) are 


Gr*( 9 )=Air‘( 9 ) 

= SjT*(9) 


for transverse current, 
for longitudinal current. 


(C5) 

(C 6 ) 
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The hclicity amplitudes can be written with the multipole amplitudes of the vector N-N* transition such as Ef± N 
Mj± N and S™ [Z± is related to the spin (Jat*) and parity (P) of N* by Jjv* = l ± 1/2 and P = (—l) i+1 ] as 




*3/2 
\NN* i 


1, 


</2 {q) = ~^ M i+ 

S^f(q) = S^(q) 


NN* 


( q ) + (l + 2)E™'( q )], 


for Jat. = / + 1/2, and 


A 


NN* 

3/2 

NN* 


(<?) = - 


- 1)(Z + 1) 


[M£**( g )+ £?£"’((/)] 


NN* i 


</2 (?) = +2 [ a + W- (?)-(i-l)C (9)]. 
S» / r(q) = S?_ N '( q ), 

for J/v* = l — 1/2. The multipole amplitudes are parametrized as 

'(A|/?-) 2 + to: 


m!± n 


(q) = 


E?± N '(q) = 


q 

m. 

A ^ 

m 


(A%™) 2 


tNN* , 


(l±l) 


) 


(A ^T') 2 + W 


M x 

2 N [2+(Z±l)/2] 


Z± (Q )> 


( A -“) 2 + g 2 


oNN* 

5 Z± 


(«)=( — 
’ TO 7 


(Z±l) 


\ZV 

(A^-) 2 + to 2 

(A-“-) 2 + g 2 


/ 


taNN* 

h l± 


(Q 2 ), 


2 N [2+(Z±l)/2] 


~iNN 

>l± 


(Q 2 ) • 


(C7) 

(C8) 

(C9) 

(CIO) 

(Cll) 

(C12) 

(C13) 

(C14) 

(C15) 


One exception for the above parametrization is applied to the first bare state in -P 33 for which we use the following 
forms: 


^3/2 P 33 (9) = -^ 3 / 2 (Q 2 )^ [G^ l (Q 2 ) - (1 - 1V)G| L (Q 2 )] , 

^/ 2 P 33 (< 7 ) = -x Al/2 (Q 2 ) l -A [G^ l (Q 2 ) - (1 +1V)G| L (Q 2 )] , 

Sl% P33 (q) = x Sl/2 {Q 2 )B G^*(Q 2 ), 


with 


A = eW\q\ 


rriA+m N \ / 1 E N (q R ) E N (-q) + m N 


1 


m N 


EN{~q) ( mA+m N ) 2 + Q 2 


iV = 


4VE|q| 2 


B = -e 


E N (-q) + mN (mi - mjv) 2 + Q 2 
1 1 Pzv(te) 41 T|z7| 4 


/ to a + m N 

2qR mN ^2 En{— q)[EN{—q) + tun] V 2toat 
1 


[(toa + toa) 2 + Q 2 ][(m A - miv) 2 + Q 2 ] 

G ? (Q 2 ) = G € (0) ( (! + a Q 2 ) exp(-6Q 2 ), 


l + Q 2 /0.71(GeV/c) 2 


(C16) 

(C17) 

(C18) 

(C19) 

(C20) 

(C21) 

(C22) 


with £ = M, E, C and Gf^(O) = 1.85, G| L (0) = 0.025, G^ L (0) = -0.238 [IM Hog; a = 0.154 and b = 0.166 
(GeV/c) -2 Il03f . The cutoff A|/?- and the coupling constants M£[ N ( Q 2 ), E[± N ( Q 2 ) and Sf±* ( Q 2 ) in Eqs. (IC13I) - 
(IC15I) . as well as the factors a:A 3 / 2 (Q 2 ), xa 1 / 2 (Q 2 ) and xs 1 / 2 (Q 2 ) in Eqs. (IC16MC18I) are determined by fitting data. 
In Ref. [32|, we have done a combined analysis of nN, *yp —> nN , pN, KA, KT, reaction data, and fixed these parameters 
at Q 2 = 0 for the proton target. The numerical values for these parameters are also presented in the reference. 
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TABLE II. Parameters defined in Eq. m for the proton target. The first column shows bare N* specified by their partial 
waves; (1) [(2)] indicates the first [second] bare state in the specified partial wave. For F^ N t(Q 2 )=M^ 1 ( Q 2 ), E[± ( Q 2 ), and 
§i± ( Q 2 ), the unit for a parameter is [1CP 3 4 5 6 GeV -1 ^ 2 (GeV/c) -2n ]. Only for P 33 (l) for which we present F^ N »{Q 2 )= 
xa 3/2 {Q 2 ), xa 1/2 (Q 2 ), and xs 1/2 (Q 2 ) from the left to right, the unit for a parameter c£ is [(CeV/c) _2n ]. 


F V 

pN* 

„P 

c 0 

„P 

C 1 

C 2 

(Q 2 ) 

„P 

C 3 

_P 

C 4 

C P 

C 5 

„P 

c 0 

c p 

C 1 

e,N* 

E l± 

c p 

C 2 

(Q 2 ) 

cP 

c p 

c 4 

•S 

c p 

c 0 

„P 

C 1 

zN* 

S l± 

c p 

C 2 

(Q 2 ) 

„P 

C 3 

c p 

C 4 

c p 

C 5 

Sll(l) 







-47. 

-429. 

754. 

-648. 

237. 

-31. 

-0.5 

9.0 

-22.3 

21.5 

-8.3 

1.1 

Sn(2) 

— 

— 

— 

— 

— 

— 

-13. 

113. 

-182. 

141. 

-50. 

7. 

0.29 

-2.55 

-0.74 

3.64 

-1.91 

0.29 

J=li(l) 

4. 

-101. 

60. 

-7. 

- 1 . 

— 

— 

— 

— 

— 

— 

— 

-0.7 

6.6 

-16.0 

13.0 

-4.6 

0.6 

p ll(2) 

66. 

-303. 

530. 

-417. 

143. 

-18. 

— 

— 

— 

— 

— 

— 

-4. 

-455. 

307. 

-86. 

8. 

— 

p ll(3) 

-28. 

161. 

-873. 

971. 

-420. 

62. 

— 

— 

— 

— 

— 

— 

13.4 

-32.0 

35.4 

-18.8 

4.7 

-0.4 

P 1 3 U) 

29. 

-133. 

-498. 

629. 

-240. 

30. 

-11. 

-158. 

437. 

-385. 

136. 

-17. 

1.7 

7.5 

11.2 

-23.0 

11.0 

-1.6 

p 13(2) 

-57. 

339. 

-689. 

544. 

-190. 

24. 

-101. 

802. 

-2074. 

1799. 

-638. 

80. 

-24. 

203. 

-182. 

59. 

-6. 

— 

Di 3 (1) 

20. 

-38. 

205. 

-176. 

59. 

-7. 

44. 

-503. 

1526. 

-1253. 

414. 

-49. 

36. 

28. 

-246. 

429. 

-210. 

31. 

Di 3 (2) 

-2.8 

15.3 

-45.9 

42.0 

-15.1 

1.9 

-7. 

-195. 

420. 

-348. 

121. 

-15. 

40. 

-103. 

131. 

-59. 

9. 

— 

Di 5 (1) 

-2.7 

15.6 

-3.7 

9.0 

-5.5 

0.9 

2.8 

-18.0 

33.1 

-24.9 

7.8 

-0.9 

19.4 

-65.2 

98.6 

-67.6 

20.8 

-2.4 

Di 5 (2) 

-19. 

-365. 

920. 

-803. 

294. 

-38. 

-4. 

-34. 

103. 

-94. 

35. 

-5. 

3.2 

-24.3 

53.9 

-48.0 

17.1 

-2.1 

p 15(l) 

1.56 

-5.27 

7.86 

-5.42 

1.71 

-0.20 

1.6 

0.7 

17.0 

-24.9 

11.1 

-1.6 

4.0 

-13.5 

21.8 

-10.6 

1.6 

— 

p 17(l) 

0.75 

-4.30 

5.53 

-2.76 

0.44 

— 

-0.13 

0.67 

-0.73 

0.26 

-0.03 

— 

1.22 

-0.24 

-1.74 

1.73 

-0.61 

0.07 

Gl 7 (l) 

3.0 

-27.6 

32.2 

-13.5 

1.9 

— 

3.6 

-16.7 

20.7 

-9.5 

1.4 

— 

12.2 

-32.9 

34.3 

-14.1 

2.0 

— 

Gi 9 (1) 

-0.3 

-21.4 

17.5 

-7.1 

1.1 

— 

-0.11 

2.78 

-3.60 

1.70 

-0.25 

— 

0.72 

-1.73 

3.58 

-1.95 

0.32 

— 

Hist 1 ) 

0.06 

— 

— 

— 

— 

— 

-0.06 

— 

— 

— 

— 

— 

— 

— 

— 

— 

— 

— 

S 3 l(l) 

— 

— 

— 

— 

— 

— 

260. 

-1903. 

4324. 

-3470. 

1189. 

-146. 

-0.61 

2.67 

-3.59 

8.07 

-4.36 

0.69 

S 3 l(2) 

— 

— 

— 

— 

— 

— 

-256. 

-6013. 

17817. 

-17166. 

6567. 

-872. 

271. 

-181. 

-88. 

175. 

-78. 

11. 

p 3l(l) 

22. 

336. 

-673. 

499. 

-161. 

19. 

— 

— 

— 

— 

— 

— 

-2.12 

-2.01 

0.79 

-0.13 

— 

— 

p 3l(2) 

-73. 

-11741. 

30223. 

-27142. 

9778. 

-1236. 

— 

— 

— 

— 

— 

— 

97. 

-393. 

754. 

-560. 

180. 

-21. 

P 33 W 

0.99 

0.54 

-0.40 

0.16 

-0.03 

— 

1.22 

1.43 

-3.62 

3.51 

-1.43 

0.20 

-0.2 

3.5 

-12.9 

13.1 

-5.0 

0.7 

p 33i 2 ) 

94. 

-352. 

576. 

-435. 

149. 

-19. 

-1.22 

-1.58 

6.42 

-5.28 

1.79 

-0.22 

-0.25 

2.98 

-1.63 

0.27 

— 

— 

^ssC 1 ) 

-3.3 

-28.3 

69.3 

-55.1 

18.3 

-2.2 

19.9 

29.7 

-74.1 

67.2 

-28.3 

4.3 

0.21 

-1.60 

3.51 

-3.06 

1.10 

-0.14 

D 33 (2) 

0. 

-192. 

486. 

-448. 

167. 

-22. 

109. 

-93. 

12. 

178. 

-100. 

15. 

184. 

-108. 

-246. 

433. 

-191. 

26. 

d 35 (i) 

-3.13 

-1.06 

4.01 

-0.62 

-0.09 

— 

-0.59 

0.39 

0.31 

-0.26 

0.05 

— 

4.83 

-6.69 

0.16 

3.56 

-1.77 

0.25 

n 35 ( 2 ) 

-10.3 

-1.6 

-2.3 

4.0 

-1.0 

— 

-9.01 

8.03 

-8.87 

4.70 

-0.80 

— 

39. 

143. 

-44. 

1 . 

1. 

— 

p 35(l) 

-0.44 

1.41 

-1.29 

0.49 

-0.06 

— 

-0.73 

0.52 

-0.07 

-0.25 

0.07 

— 

5.6 

6.0 

-10.9 

4.7 

-0.6 

— 

p 3 5(2) 

-0.5 

-26.1 

42.1 

-21.2 

3.4 

— 

40.1 

-76.4 

69.4 

-27.3 

3.8 

— 

0.02 

0.11 

-0.43 

0.16 

-0.02 

— 

p 3 7(l) 

0.48 

-2.84 

4.59 

-3.15 

0.98 

-0.11 

0.02 

0.06 

-0.10 

0.05 

-0.01 

— 

-0.04 

0.20 

-0.27 

0.13 

-0.02 

— 

^37 ( 2 ) 

11.2 

34.7 

-16.2 

12.5 

-3.2 

— 

-1.8 

-57.3 

83.5 

-43.6 

7.2 

— 

6.2 

-17.8 

26.6 

-12.2 

1.8 

— 

G 37 (l) 

-0.30 

0.63 

-0.58 

0.42 

-0.09 

— 

0.89 

-2.53 

1.95 

-0.65 

0.08 

— 

3.1 

-7.2 

11.9 

-5.2 

0.7 

— 

g 39 (i) 

0.2 

-9.2 

10.2 

-4.2 

0.6 

— 

-0.03 

-0.37 

-0.45 

0.52 

-0.12 

— 

0.89 

-1.48 

1.78 

-1.24 

0.40 

-0.05 

n 39 (i) 

-0.09 

1.44 

-2.16 

1.07 

-0.17 

— 

-0.22 

-1.44 

2.16 

-1.07 

0.17 

— 

— 

— 

— 

— 

— 

— 


TABLE III. Parameters c ™ defined in Eq. (1431) for the neutron target. The parameters for 1= 1/2 N* are presented while those 
for 7=3/2 N* are the same as those for the proton target. The other features are the same as those in TABLE HT1 


F V 

F nN* 

r n 

C 0 

C 1 

C 2 

(Q 2 ) 

c 3 

c 4 

c 5 

c 0 

C 1 

j?N* 

E l± 

c 2 

(Q 2 ) 

n 

c 3 

c 4 


c 0 

C 1 

oiV* 

S l± 

c 2 

(Q 2 ) 

c 3 

c 4 

c 5 

Sll(l) 







28. 

730. 

-1439. 

1092. 

-344. 

39. 

180. 

-944. 

1672. 

-1067. 

279. 

-25. 

Sn(2) 

— 

— 

— 

— 

— 

— 

38. 

-250. 

456. 

-343. 

114. 

-14. 

16. 

-67. 

105. 

-69. 

20. 

-2. 

p ll(l) 

-30.2 

25.3 

42.9 

-59.2 

23.0 

-2.9 

— 

— 

— 

— 

— 

— 

-13. 

156. 

-282. 

188. 

-47. 

4. 

p ll(2) 

-118. 

222. 

-345. 

256. 

-83. 

10. 

— 

— 

— 

— 

— 

— 

-11.9 

34.5 

-42.4 

23.1 

-5.7 

0.5 

p ll(3) 

107. 

-527. 

939. 

-701. 

230. 

-27. 

— 

— 

— 

— 

— 

— 

2.4 

-32.7 

-1.1 

27.8 

-13.8 

1.9 

P 1 3 U) 

-2. 

-76. 

126. 

-98. 

36. 

-5. 

5. 

-85. 

115. 

-78. 

27. 

-4. 

-13.2 

33.8 

-32.4 

13.9 

-2.9 

0.2 

p 13(2) 

124. 

152. 

-30. 

-26. 

7. 

- 0 . 

34. 

70. 

-186. 

116. 

-21. 

- 0 . 

25.4 

-20.7 

-21.3 

39.6 

-16.9 

2.3 

Dl 3 (l) 

-35. 

216. 

-521. 

392. 

-120. 

13. 

-21. 

343. 

-731. 

564. 

-180. 

21. 

139. 

-379. 

426. 

-113. 

-3. 

3. 

Dl 3 (2) 

20. 

-109. 

193. 

-139. 

44. 

-5. 

-14. 

90. 

-164. 

119. 

-38. 

4. 

6.6 

-48.5 

70.6 

-42.6 

11.6 

-1.2 

Di 5 (1) 

22. 

-117. 

208. 

-154. 

50. 

-6. 

0.9 

-8.0 

17.4 

-13.3 

4.3 

-0.5 

6.1 

-24.6 

35.3 

-22.9 

6.9 

-0.8 

Di 5 (2) 

-30. 

154. 

-271. 

193. 

-60. 

7. 

-2.1 

11.3 

-20.0 

15.1 

-5.0 

0.6 

13.8 

-49.4 

67.6 

-42.8 

12.6 

-1.4 

p 15(l) 

-1.8 

14.6 

-18.5 

10.0 

-2.7 

0.3 

0.6 

-10.6 

44.9 

-44.0 

16.2 

-2.0 

40.5 

-88.5 

48.1 

15.3 

-14.9 

2.5 

^17( 1 ) 

-1.37 

1.45 

0.80 

-4.85 

2.84 

-0.46 

3.9 

-16.3 

26.9 

-17.3 

4.9 

-0.5 

-2.0 

-5.3 

16.0 

-14.0 

4.9 

-0.6 

g 17( 1 ) 

2.77 

-8.29 

7.62 

-6.36 

3.11 

-0.52 

-19. 

112. 

-202. 

147. 

-47. 

6. 

-4.32 

2.07 

-5.28 

5.77 

-2.16 

0.26 

^ioC 1 ) 

-2.7 

14.4 

-25.1 

18.0 

-5.7 

0.7 

0.13 

0.73 

-2.75 

2.06 

-0.54 

0.04 

16.8 

-46.8 

49.0 

-23.7 

6.0 

-0.6 

H 1 9 (1) 

0.84 

-3.22 

4.25 

-3.14 

1.17 

-0.16 

-3.3 

20.0 

-36.5 

26.5 

-8.3 

1.0 

10.7 

-28.9 

29.6 

-13.7 

3.4 

-0.4 


Appendix D: Parameters in Eq. (1431) 
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TABLE IV. Parameters defined in Eqs. (IC13I) - (ICT5I) for the neutron target. The parameters for 1= 1/2 N* are presented 
while those for 7=3/2 N* are the same as those for the proton target and are given in Ref. ?32l. 
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